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functions of the National Bureau of Standards are set forth in the Act of 
Congress, March 3, 1901, as amended. These include the development and 
national standards of measurement and the provision of 
means and methods for making measurements consistent with these standards: 
the determination of physical y= Bere = and properties of materials; the develop- 
instruments for testing materials, devices, and structures; 
advisory services to government agencies on scientific and wry problems; 
invention and development of devices to serve special needs of the Government: 
and the development of standard practices, codes, and specifications. The 
work includes basic and applied research, development, engineering instru. 
mentation, testing, evaluation, calibration services, and various consultation 
and information services. The Bureau also serves as the Federal technica! 
research center in a number of specialized fields. The scope of activities of the 
National Bureau of Standards is —— in the following listing of the divisions 
and sections engaged in technical work. 
Washington, D.C. 
Resistance and Reactance. ~~) Electrical Instru- 
Measurements. Dielectri 


Photometry and Colorimetry.  Retract ctometry. Photographic Re- 
nes Engineering Metrology. Mass and Scale. Volumetry and 
mg Temperature Physics. Heat Measurement. Cryogenic Physics. Equa- 
tion of State. Statistical Physics. 

Radiation Physics. X-ray. Radioactivity. Radiation Theory. High Energy Ra- 
diation. Radiological Equipment. Nucleonic Instrumentation. Neutron Physics. 
Analytical and Inorganic Chemistry. Pure Substances. S emistry. 
Solution Chemistry. Standard Reference Materials. Applied Analytical Research. 
Sound. Pressure and Vacuum. Fluid Mechanics. Engineering 
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‘Mechanics. Rheology. Combustion Controls. 


Organic and Fibrous Materials. Rubber. Textiles. Paper. Leather. 
Testing and Specifications. Polymer Structure. Plastics. Dental Research. 
Metallurgy. Thermal Metallurgy. Chemical Metallurgy. Mechanical 
Metallurgy. Corrosion. Metal Physics. Electrodeposition. 
Mineral Products. Engineering Ceramics. Glass. Refractories. Enameled 
Metals. Crystal Growth. Physical Properties. Constitution and Microstructure. 

Research. Structural Engineering. Fire Research. Mechanica! 
Systems. Organic Building Materials. Codes and Safety Standards. Heat 
Transfer. Inorganic Building Materials. 
Applied Mathematics. Numerical Analysis. Computation. Statistical Engi- 
= Mathematical “9g 

C En apart Group. Components and 


os on Radiometry. Solid State Ph Electro 
Physics. Atomic Physics. * aay 4 


Instrumen Electronics. Electron Devices. Elev'tronic 
Instrumentation. Mechanical Instruments. sess Instrumentation. 
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Molecular Spectroscopy. Molecular toe Mass Spectrometry. Mo- 
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e Office of Weights and Measures. 
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rain Effects. Radio-Meteorology. Lower Atmosphere Physics. 


Radio Standards. High Frequency Electrical Standards. Radio Broadcast 
Service. Radio and Microwave Materials. Atomic F: ency and Time Interval 
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wave Circuit Standards. 
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Optimal Approximation for Functions Prescribed 
at Equally Spaced Points 


H. F. Weinberger 


1. Introduction 


Many probl mis in numerical anolvsis Cul be re 
duced to approximating the value F(z) of 
linear functional / operating on an unknow! 
The 


set of data concerning wu 


vIVeH 
element 
vu of a linear vector is to 


space approxine tion 


terms of a finite 
values F)(u om 1) of N linear 
Kor example, 


be made ith 
Thus, 


functionals a 


thre 


cting on wmavy be give 


the F(a) may be values of the function w at certain 
points If Fu is the value of w at anothet pot 
£, we have the problem of linear interpolation If 
F(a) is an intezral of 4, We have the problen olnu 


merical quadratures If F(a) is the value of a deriv 
ative, we have numerical differentiation 

| was shown by \l (Lolomb and the author that 
in order to obtain a finite interval in which the value 
fu) must le, one must be rive the value of at least 


The 
Is given a bound ior a 


(CLolomb 


one nonl near tune tional operating On UW Sill 


plest ease is that in which one 


quadrath functional iu In this case 


and the nuthor showed how to obtain the exact 
interval in which Fw) must le whe the values 
ki(u / 4), and ‘u) are given That is, 


which are attained 


bounds lol / iM 


sutisiving the 


uppel ind lowel 


for some elements / viveln conditions 
ure found The 
quires the inversion of a matrix depending 


functionals F’, / / 


construction of these bounds re 
ipon the 


and Mou 


In this paper we restrict our attention to a very 
simple cause We deal with a function u(r) of a 
single variable on the interval 0.1) The given 
functionals /, are the values wu V) of wat the A ] 
equally spa ed pomts \ U, \ The 


Pur tional is taken to be the integral of the 
square of the &™ derivative of u 

We think of the number of pomts V as large, 
the number & of the derivative will usually be 
say two or three The 
size .\ $y making use of the equal spacing of our 


quadrati 


while 
small 
matrix to be inverted is of 


10. 196] 


pots wi shal! reduce its inversion to that of a 
j | j 1) matrix Thus the 


st formulas for interpolat ol 


ob- 


problem ol 


taining be quadratures, 


and numerical differentiation Is ice inanageable 
even whel the number ol points nvolved becomes 
leat re 


When / f ide. our results vield us special 


cases the best quadrature formulas of Sard for 
j .N<6 
2. Approximation Problem 
Let e values of the unknown function ula be 
LV ¢ il thie \ | evenly spaced points 4 \, 
) \ Let \/ be a rivel bound for the 


eral of the &th derivative of u 


We assume that \ > DI j 2 
dur problet Is tO approximat thie 


value F(u) of a 


certain linear functional - applied to u According 


to the theory in footnote 2, this is possible if and 
only if the functional F’ is bounded in the norm (2.1 
for functions vanishing at the points 7/.\ That is, 
we must assume that there is a constant ec such that 


*1 
F(v)*<e | ly 2.2) 


es differentiable functions such that 


for all & ti 


(..)=0, 0 TN 2.3 


Any linear combination ol values or 
integrals of and its derivatives up to order k-1 will 
this condition 


pointwise 


satisfy 














To obtain a best estimate for F(w) we construct 
an auxiliary function @ defined by the properties 


u™* (x) —0, G6<a<1, Nzst0,1,..., / N (2.4) 
and 
_ 1 1 
(x) u(x) ?=0,1,..., 1 \ 
u'4(0)=U!"(1)=0 f=-—,..., 2—2 (2.5) 


The function wu and its first 24-2 derivatives are 
continuous, while @**~" is allowed to have jump 
discontinuities at the points 7/N. 


The Green’s function G(x; &) is defined by the 
properties 
Oo" , ‘ oa 
— =(), 0<2<1, Nr¥N¢,0,1,...N, 
Ox 
G ( 5 E) 0 baeO,1, .. x» J \ 
\ (2.6) 
0'G 
7=0 at z=0,1,l=k, ... 2%—2 
or 
oO” ‘GE +-()° £) Oo” 'G(é 0: &) . 
" { 
Or” l Or” } J 





Again G and its first 24-2 derivatives are continuous, 


while the (2k—1)st derivative may have discon- 
tinuities at 7/N. By integration by parts we find 
that 
N : 
u(é) >> gilé)u (55), (2.7) 
& iv 
where 
g(é)=(—1)**' 
ok~1/7 i * o~177 t . 
a*-'G ( yt; 8 ) ora ( yO 8 ) 
: - - (2.8) 
on* or 
Also by integration by parts we have 
on : one OG (x; &) 
u(t)—U(¢) 4 1 — yl 22" dzx. (2.9) 
s s - or 
Applying the functional / and using Schwarz’s 


inequality, we find that 


; _ be —_— "| POG (ae) 7]? 
F(u)— F(a)? < yu) — yl “ds | Ff =“ ‘| da 
0 0 | r 


(The symbol Fn [G@ (£; n)] means that the functional 
F is applied to G@ considered as a function of » for 


u™ \*drF:{ F(G(E; n)| 


(2.10) 
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fixed & I FIG (&; n)] 
to the function F/G). 


means that Fis then applied 
We have used the property 


"T YG (E: x) OG (n: 2) 


dx, 2.11) 
Or Or 


G(é; n) 


which follows from integration by parts 
Another integration by parts shows that 


TAL, A TAL » » 
u u u dr=0 2.12) 


Hence, we can rewrite (2.10) as 


F(u)— F (ii) \?< {ar ( mine \ F\G(t:n 


2.13 


Once u and G are found, this inequality provides 
upper and lower bounds for /(w). These bounds are 
sharp in the sense that we can construct functions w 
satisfving (2.1) and having the given values u(7/N 
for which the bounds for I “w) are attained 

We write 7% in the form 


\ 


u(x) >, aA*| Nzr—i| 2.14 


The centered difference operator A* 1s defined by 


. { 2 , 
se= D> (.,))(—* te 2.15 


It is easily seen that A* No } vanishes for 
Na—i >k. Hence the sum in 
2k—-1 nonzero terms 

Clearly the function (2.14 
the required continuity properties 


has at most 


2.14 


satisiies 2.4 and has 


The coefficients a, are to be determined by the 
conditons (2.5 Thus we must have 
N+k—1 j 
S* aA* j—i|* u ( vy) j=0,1,..., N 


In order to apply the last line of (2.5), we first use 
partial summation to write 


ua 1)* » > L¥Fa, TD Na—i* : 2.17) 


where we have put a 0 for 7< k i>N-+k, 
where Ly is the Ath forw ard difference operatol 


iT 


Since A*| Nx—i|”* 0 for |Na—i|>k, (2.14) and 
(2.17) are independent of the values of a, for 7: k, 


i>N+kwhen 0<2<1. We note that D*|\Nr—i 





Is a polvnomial of degree f 1 for Nx <ior Nr>i+h 


Hence the last line of (2.5) will be satisfied if 


La 0). H j B 
1, N—k Be ae , A l 2.19 


The boundary value problem 2.5 is thus replaced 
by the svstem (2.16), (2.19) of N-+2i 1 linear 
equations in the .V+2k—1 unknowns a 


We turn to this problem of matrix inversion, 
Re mark We have assumed that ke 2 The cause 
h 1 is easily treated Let %@ be the broken linear 
function coinciding with 7@ at the points 7 V and 
having its breaks at these points Let G r:t be 
the broken 
two neighboring points of the form i/N which 


inear function with breaks at N and the 
vanishes at the pots V, and whose derivative 


decreases by 1 at ¢ Then (2.13) gives optimal 
bounds for I: u when 4 l 


3. Matrix Inversion 


We consider the svstem of linear equations 2.16 


2.19 since A* ) depends only upon ” 
and vanishes for j 2.16) is a finite dif 
ferences equation ol ordet 2h 2 We solve il by 
means of a system of 2A 2 independent solutions 
of the homogeneous equation To final these solu 


tions, we note that for anv number 
thm} 
‘ 


e*A*| j—i|” 26 


for 0 \ 
We define the polvnomial Ith ¢ 


d 
LE e?7—] 


lt is easily SeCeECT! to he ol degree / with leading 


coefficient 2 
The i can bye venerated by the recursion 


(/ y 2 y ] 
().(y 2(y ty 1), 24 
Q3(y)=2(y*' + 11y?+ 1lly+1), 
(), y 2 y' 26y 66y° 26y | J 
The coefficients of t are the coefficients of the 


finite difference equation (2.16) It can be shown 


by induction that thes ure svimmetric in the sense 
that 


This means that the zeros of (/, occur in reciprocal 
pairs When / is odd, one of the zeros is l When 
liseven, Y can be written as a polynomial in (y+ y 
Thus, the zeros of (Jox-2 CaN be found by solving an 
equation ol degree k l and a quadratic equation 

It is shown bv induction that the coefficients of 
(J, are positive and that the zeros of (/, are real and 
negative The zeros of (/, separate thos of A 

Let 7] 7 7 0 be the zeros of 


Because of 1 and 5.2 the functions a Y, 
satisfy the homogeneous equator Ss corresponding to 


2.16 We now define 


0 <k-—2 
! 2k 5.) 
¥ — ¥ b-9 
ow «) 7 = 
lt follows from the Lagrange interpolation 
formula [4] that the two definitions of ¥, coincide for 
k+-2: k—2. Using (3.1) and (3.2), we find 
that 
\ 
SS yy pA*lj—il* 5... V.p=9, NN 3.9) 


The same equation is satisfied by ¥ plus any 
linear combination of the y We add a linear 
combination such that the new function satisfies 
2.19 since y¥ vanishes for ¢ « p k}—2 these con- 
ditions are already satisfied for 1 j l 2 


There fore we add only fun tions which also satisfy 
these conditions The Lagrange interpolation for- 


mula * shows that such linear combinations are 
furnished by the functions n,,., a=! _k, where 
—_ U] 
y S 3.10 
=t (¥—1)°Y y 


Note that 


Dn, =vW,, k+2 S11) 
We let 
I ¥ >) Cati+e (3.12) 


and determine the coefficients ¢, in such a way that 


‘J. F. Stetle nterpolation (Chelsea Pr New York, N. Y., 1950 
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DiT,,=0 for i=—1, N-—k+1,..., N—1. This 
gives ¢,=6),, and the equations 
k-1 


y Ca Ww k+a+6 D Wy k+8-p 
=1 


— bop Vw +8, B=1,...,k—1. (3.13) 


a 


We let As be the inverse of the symmetric 
(k—1) (k—1) matrix Py_,,048 80 that ° 


k—1 
>) Aas Ww—248+47 - Sey; a,y ] ak l 
8=1 
(3.14) 
Then I;, is given by 
k—1 
Dts vi a Aas Ni+a LP W k—p+8 
as=l 
Sop 7 h >) Aes Nita Wn 3 t 
k+-1,..., J N+k—1, p=n0,..,, J WV. (3.38) 


This function satisfies 


(3.16) 
and 
Dit ,.=0, i k+1,...,—1, 
N—k+1,..., N—1; pus), ..., J N. (3.17) 
The solution of (2.16), (2.19) is given by 
Sop (2). 
a 2! pU\ (3.18) 


Therefore [,, is the inverse matrix for the problem 
(2.16), (2.19). 


Our problem of matrix inversion has thus been 
reduced to the solution of polynomial equations of 
degree k—1 and two in order to find the y, and the 
inversion of the (k—1)-dimensional symmetric ma- 
trix Wy—.40+28- 

Example: We consider the case k==2. 
of Q,(y) are 


The Zeros 


Yi 2—y3, 
y2=—2+73. (3.19) 
The function y, defined by (3.10) is 
0 i< 0), 
W=<4 V3, , (3.20) 
19 Yi —Yy) > 0. 
‘The fact that vn s+e+3 is nonsingular follows from the uniqueness of u le 


fined by (2.4), (2.5), and the linear independence of the functions A?*| N 
and give,a=1,.. ., &. 


1, the matrix A, has the single element 


Moreovy er, 


Since k—1 


vx’. 





- l 
py 5 V38(y yit'), i>0, 
a 
¥ ? 
| ; 
n a5 VOY Y, °) (3.21 
Thus (3.17) gives 
fac) y; *)(l—y, y 
12(1—y,; *) 
l<i.<p l<p<N-1 
\ } | Y (1 Yi; . vi 
12(1—y,**) 
P / \ # l pP \ l 
r 
l—y, *)yi~* 
\, p= N 
12(1—y; >) 
(1 y, 2 \y 
= i—0 p=—0 
12(1—y, ?*) 
} \o (4 \»” y \ | p 0) 
L 12(1—y; ti=N+1, p=N 
») 


4. Bounds 


return to the consideration of section 2 


We how 


The function 7% is given b: 


N+k—1 
Ua y & a,A” Na , ? . (4.1 
k+1 
Ni 
with 
N : 
Dl, u( +) 
a I u _ |} +2 
\ 
Thus, we may write 
N 
u (zr) >> 9; (x) u( ) 4.33 
J { 
where 
N 1 
g, (a S r,,A*| Nx—7 1.4 
i 11 
\ I 


is the optimal approximation function corresponding 
to u(i/N)=6 


In order to find the bounds (2.13) we need the 
integral of 7? and G(z,¢ Integrating by parts 
we find that 

*] N 
i"? d2—(—1)* 3+ u( 4.) (a tan 
. Jue N 


y denotes the discontinuity in the function 
with the 


Applving 


W here 

ul IN, 
Be 0) 
we find that 


convention % 0 7} 
partial summation to (4.1 


\ 
u (a >» Vu Aa 1.6 
Consequently 
— N 
j, weds N*-1(2k—1)! SY uw ( ay) A*a 
»\ 2k—1)! 


Su(i)u(f)atry 7 


with the convention that 


AT ta fol 0.4 \ 1S 
We now construct the Cireen’s function G rs 
We begin with a function //(2,£) having the propet 


the 
the end 


jump at g and satisfying the condition that 
derivatives of orders &. 2k-2 vanish at 
points. Let p be any integer satisfying 


Fu F(A Vu u( ‘ ) 
V2—2(—1)*N 
2(2k—1)! ~ 
\ 
where we again use the convention LS If I Is 


local in the sense that /'(w) only involves the values 
of uw in the neighborhood of a point, the sums in the 
second term on the right involve at most 2A 1 values 


of ¢ and 7. exam ple : Letk—2 l,, is given by 3.22). 
We find that 
> 4—6 Noa | Na 
A’ Na 2{2—!Na 
{} 
Va i<l 
l Na 2 
Va 2 $15 


It follows from the Lagrange interpolation formula 
that the function 


see footnote 4 


Ha rt >| b 1.10 
where 
V: 2k—1 
} 7 | : ( : ) 
fp MM jf} i 
(e—\-) om (&—"Q-) (41D 
has the property 
. , j 2} | 
HT 2:3 (0) for «2 \ \ 4.12) 


Therefore, the function 


{ N 
Gila Hi(x;¢ —>) } r) HT | =F. } 
(4.13) 
has all the properties (2.6 Thus, it is the Green’s 
function We not that the sums in (4.4) and (4.13) 
involve at most 2/ 1 terms for each z and &. 
The bounds (2.13) for F(w) are now given ex- 
plicitly by 
. } 
) hb , d 
2k—1 seriu(x)u(Z) } 
. . 
SP, F(4*| Nz—i r[u( x8) ]} (4.14) 
Thus the interpolation function g,(2) is given by 
qg,(a 12 q l Va Vz QMi(q+2 Va)l 
Vzr-+-1 ql } q+ Vr)* 
Va j/) 0 | 16) 
where the integel q is defined by 
If ] NE \ l. we let Pp be the inteve! such that 
p+i1< Né&< p+2 4.18) 











Then we have 





4 () Na P; Na > po 3 
(Nz P) >. = 
NE—p ps Nespt! 
, (Nir—p 3( Nr | : 
: ___ 3N*H(z; &) — P) “Ate P 7: pt+1l<Nar< NE ™ 
(NE—p)(NE—p—1)(p+2—Ne)(pt3—Ne) ) ** = s (4 
12 "+ 2(n+9 To» 
(Pro Nz) OY é Nz) Ni Nr<p+2 
p+3—Né p+2—Ne¢ ; 
(p+3—Nz) a ti ; 
. _" Pres Na Pp D- 
. p To N¢ 
If 0< NE< 1, the definition (4.18) gives p 1, so 
that the function 17 defined by (4.19) does not vanish 
r=0. his case » si subtract the | , : ; : 
ra ts a yee we apy swtire | from the right-hand side of (4.19) with p=N—2 
unction poo 
Thus we can evaluate the bounds 4.14 explicitly. 
In the special case of linear interpolation we have 
(NE+1)—'A*| Nz—1)8 (4.20) 
9 | 
é F(a) u(t). $1.22 
from the right-hand side of (4.19) with p | to 
obtain H(z, &). Similarly, if N—1< NECN, we If 
subtract 
l<¢: Ne q+: N—1, (4.23 
pe ‘ ; . 
(N+ 1—Ne)-'A*| Na—(N+1) (4.21) | . 
2 | we find 
re N+1 N ; j 
FAF(H(2;§))|— >> >) 1, F(A*| Nz—i\’)F [7 ys | 
t 1 j=0 d 
N-*(Ne—g)*(q+1 vor{ i—s7 a 3) | (a 2—Ney*l—yyy—y 
+ 2(q+2—N¢)(No&—q4+1)—y; (l—y;?% yyy! 
(Ne q- 1)?(1 Yi” (1 Yi i |} (4.94 
where bound for the square integral gives considerably 
oF 2—3. less information than a bound for the maximum 


The first factor on the right of (4.14) does not 
approach zero as N—>@ unless \/* happens to be the 
exact value of the integral of u“”. Thus the differ- 
ence between the best upper and lower bounds in 
the linear interpolation problem with N+ 1 equally 

7 
spaced points and with a given bound for u’ 
e 0 
is of the order N If a uniform bound for |x’ 
is given, one can obtain bounds for u(¢) which 
differ by a term of order N~*. This shows that a 


"dx 


The problem of finding best bounds when the maxi- 
mum of ju“?! is bounded is much more difficult than 
the problem treated here. 


The author wishes to thank W 
a multitude of suggestions in connection with this 
paper. 
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Truncations in the Method of Intermediate 
Problems for Lower Bounds to Eigenvalues 


Norman W. Bazley and David W. Fox' 


Februar 


10, 1961 


Iw ew procedures are developed for determining lower bounds to the eigenvalues of 

ear operators The method ire based on he theory of semibounded self-adjoi 
ito n separable Hilbert space Computation of the ower bounds is reduced to the solutior 
f matrix probiems The procedures nave immediats application in the estimation rT 
genvalues and eigenvectors of differential operators appearing in quantum mechanics 


1. Introduction 


this we two 
to eigenvalues of 
bounded self adjoint operators These procedures, 
while phrased in the of Hilbert 
result in useful computational methods for the esti- 


mation of eigenvalues and eigenvectors of differential 


rive 
bounds 


In papel new procedures for 


determining lowe semi 


language space, 


operators occurring mm Classical and quantum me 
chanics Our procedures ure based on the method 
of intermediate problems 

Intermediate problems as introduced by A. Wein- 
stem Z| rive lowe! bounds by changing the bound- 
ary conditions of differential operators By relaxing 
the a solvable base 
proble m that vives rough lower bounds to the eigen 
values He then links the base problem to the viven 
problem by a of intermediate problems 
which can be terms of the base problem 
and which in prove the lower bounds 

Years H. Weyl [8] 
lower bounds can also be obtained by changing the 


This fact by N. Aronszajn [1] 


to construct intermediate problen s for those cases 


boundary conditions he obtains 


seq ue nce 
solved in 


ago pointed out that rough 


operatlol was used 
in whi h i base pr ble n can be obtan ed by neglect- 
a positive term in the operatlol In the same 
papel (ronszayn svstematized the method of inter 
mediate problems for completely continuous oper- 
ators in Hilbert As Aronszajn pomts out, 
the solution of the intermediate problems in his 
formulation requires the determination of the poles 
and zeros of a meromorphic function usually given 
as an infinite sum 

Recently H Kk Weinberger that it 
to overcome the principal difficulties that 


neg 


space 


Is 


showed 


possible 


arise in the intermediate problems ol Weinsteint 
Weinberger’s result, which is a special case of his 
more reneral theory 6}, fives lowe! bounds to the 


eigenvalues by the computation of the zeros of a 
polynomial 
In this paper we first sketch Aronszajn’s method 


of intermediate problen s for changes in the operator 
Qur presentation allows for the presence of continu 


\ ratory, The Johns Hopkins University r Sprit 


td 


I ket 


spectra and is in a form suitable for application 
many quantum mechanical problems In the 
second section we show how it is possible to modify 
the method so that lower bounds can be obtained by 
elementary calculations. These results give formu- 
las parallel to some of those of Weinberger 16] A 
used to 


ou 


to 


feature of our procedure is that it can be 
calculate lowe! bounds to the energy levels of most 
atomic systems. The last extends another 
method of one of the authors 

Part of the work presente d here was s arted while 
the authors were in the Institute for Fluid Dynamics 
and Applied Mathematics, University of Marvland, 
and was summarized in [3] 


section 
>] 


2. Method of Intermediate Problems 


H isa separable Hilbert space 
with the inner product and that A is a if- 
adjoint operator with domain Y, in #” We will 
A is bounded below and that the first part 


Let us suppose that 


sc 


7 


USSUTLLE 


of its spectrum is discrete. The eigenvalues and 
eigenvectors in this part of the spectrum will be 
denoted by A, and u, l 4 respectively ; 


the first limit point of the spectrum will be called 
\ Thus we have 


Au,—dr,u,=0, pax] 2. 1) 

and 
Ai S Ao‘ = (2) 
Suppose further that A can be written as the sum 


of two operators, 


where A® is a self-adjoint operator with a known 
spectral resolution KK and A’ is symmetric and 
positive The domains of A and A’ ure denoted by 
y*, and 4 as respectively, so that Y, 44 44 


The positivity of A’ is expressed by 


>{) (4) 


(A’u.u 
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for every u in the domain of A’. Finally we assume 
that the first part of the spectrum of A° is discrete, 
and, listed in increasing order, the eigenvalues are 


MSAS.... 2) 


The corresponding eigenfunctions are u?, wu, 
Since A® <A, it follows® [4, p. 214] that 


nN <A,. (v=1,2,... 
and 
r°, <i... (6) 


These inequalities go back to H. Weyl [s]. Thus 
the known eigenvalues of A® give rough lower bounds 
for the initial eigenvalues of A. For this reason the 
operator A® is called a base operator and the eigen- 

value problem for A® is called a base problem. 
The procedure of intermediate problems links the 
base operator A° to the given operator A by a sequence 
A*} of intermediate operators. These operators sat- 

isfy the inequalities 

A®°s AS A*t' <A, (k 


(7) 


bo 


and consequently their eigenvalues give intermediate 
lower bounds to those of A. 

In order to define these intermediate operators, 
we introduce an auxiliary Hilbert space #7’. .#” is 
just the completion of Gy in the norm generated 
by the inner product [u,v] defined by 


[u,v] =(A’u,r). (S) 


For simplicity we will assume that <A’ is positive 
definite; the indefinite case requires only a slight 
extension of the theory. Let {p,} be a linearly inde- 
pendent set of vectors in %,4; note that since A’ 
is positive definite the sequence { A’p,} is also linearly 
independent. Define P* to be the projection in .# 
on the span .#* of the first k elements of {p,}. For an 
arbitrary vector v in 4%’, P*v is given by 


, 


—— 
t, 


P*r= S~* [p, pilb Ps, (9) 
l 

where the constants 6,; are the elements of the 
matrix inverse to that with the elements (A’p,, p;). 
Since the subspaces .* are increasing with /; 1.e., 
2 pate a (F=1,2,...), (10) 
the corresponding projections satisfy the inequali- 

ties, 


0<P*<PHi<] (11) 


; 


where J is the identity operator in #’. We now 
form the operator A’P*, which is well defined on 


/ ¢ 
1°<A means 74°=> 4 .and (A%u,n)< 1u,u) for every in 7 4 


Y,. For an arbitrary vector win Gy, A’P*wis given 


by A’P*w= >> (w,A’p,)b,,A'p,, (k=1, 2,.. .). 2 


In the form given on the right in eq (12) A’P* has 
an obvious extension to all of 4. and there itis a 
symmetric completely continuous operator. Hence 
forward we will regard A’P* as thus extended 

It follows from eq (11) that in _# the operators 
A’P" satisfy the inequalities, 


0<A’P*<A’P* <A’, k=1, 2, 13 
The intermediate operators A*® are now defined by 


A®'=A°+A’'F", pant, 2 14 
k is called the order of the intermediate operator 
since the operators A’P’ are sVimimetric and bounded, 
it follows [5, p 301] that the operators A* are self 
adjoint and have the Same domain as A And 
moreover, since A’P’ is for each h a completely 
continuous operator, each of the operators A* has 
exactly the same limit points in its spectrum as does 
A [5, p. 367]. Krom eq (13) it follows directly that 
AS AS ATCA, (@==1, 2: 
These inequalities ure just those stated in (7), and 
they show that the operators A* defined by eq (14 
are in fact intermediate operators This means that 
in the initial discrete part of the spectra the eigen- 
values satisfy the inequalities 


N<AE<EHISA,, (v=1,2,...;@), 15 


where the superscripts on the eigenvalues indicate 
the order of the intermediate operator to which 
they belong. 

We now turn to the problem of determining the 
eigenvalues and eigenvectors of the intermediate 
operators <A’ We will confine our attention to 
those eigenvalues which lie below the first limit point 
dX”, common to all operators A*. 

The eigenvalue equation for A“, 


A*u Au 0), 16 


can be written in a more useful form by using the 
defining eq (14) for A*. This gives 


A°u—ru A’P*u. 17) 


and by using eq (12) we obtain 


A°u—du >) (u,A’p,)b,)A’p,. (18 
l 


For convenience we designate the quantities pe 
so that the equations to be satis- 


(u,A’p \h by ay, 
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fied become These together with eq 24 determine the constants 
a The necessary and sufficient condition for the 
Au hu A’ p,, (19 existence of nontrivial solutions is the vanishing of 

the determinant of the coefficients; i.e 


and /Q R A’p,,A’p A'p 


det 0 97 ) 
‘a j A’p u,A’p . / > ~ ° 20) , , ) 


S 
ya \’p 
Let us assume that (27) holds; that is, that A, 1s an 
lf fis to be an eigenvectol corresponding to a eivenvalue ol A us well is ol A The eigenvectors 
value \ which is not in the spectrum ol ig then u of A corresponding to A, are then determined in the 
may be written following way: The eqs (24) and (26) have linearly 
indepe ndent solutions a>. @ a oO 
u > @ RA’ p 2 Pe 4 n), where nis the nullity of the coefficient 
matrix displayed in (27 sv (25) these give rise 
where R is the resolve! | operatlol lol l On VO tors . that are eigenvectors ol A for the 
ah ) value A lhese are also linearly independent, for if 
Th eqs 20) become 
( are constants such at & Cet 0, then bv the 
expression 25) for uv we must have 
Sa RA’ p,,A'p)=0 2 I 22 
Re S Daz A’p, +O. S a 0. 
All of the vs cannot vanish if wis to by nontrivial 
' , Is) 
therefore \ must be such that the determinental S 
y ut n ) 
CQUAtLO But since Jt.» is a one-to-one operator on the orthog- 
p ie 1’ p,, 1’ p () aa onal complement ot. 7 , eg 2S) 1 phic s the separate 
: Statements 
Is satisfied. 
Kquations 2 have solutions for each value \ that . 4 , 
Is not in the spectrum of A° and that satisfies eq (23 2 ( 20° a ).A’p 0 >> ( 26 a ) 0 


The number ol solutions corresponding to ene fh such 
24) 


root is just the nullity of the coefficient matrix and, 


according to eq (21 these solutions give the eigen The linear independence of the vectors A’p, and of 

vectors ol \ ussociated with that root These the vectors imply that 

eigenvectors are cle rly linearly independent since ? 

R Is one-to-one oO! Hf 

Ion euct value \ that is an eigenvalue ol thre > . 0) 1,2 | / 

operatol A’, say \ the procedure Is somewhat 

lithe rent Let y be the charac Leristi subspace But sine the vectors at. a — are 

of A° at the eigenvalue ) Suppose that it Is | Jinearl independent, each of the constants C, must 

spanned by the eigenvectors . Its | vanish, and thus the vectors wv’ are linearly inde- 

clear that if w is an eigenvector of A corresponding pendent 

to A then in addition to the eq (19 and (20) the As wi observed previously, the part of the 

further restrictions spectrum ol A* that lies below \ is discrete The 
eigenvalues there are completely determined from 

Sal A’p 0. ] 1 | ) ] (24 eqs 23) and (27). and the eigenvectors are found 

—, : from corresponding solutions of the linear equations 
These eigenvalues vive intermediate lowe! bounds 

Hiust be satisfied lf we assume that eq 24 nre to those in the first part of the spectrum of A 

satisfied “wus Given hy For the suke ol completeness and late use, Wwe 
sket« hy il proof of convergence of the eigenvalues of 

\ » A'p X 25, the intermediate problems The conditions that we 
1 =k use are those stated by Aronszajn i1] We assume 


that A and A have completely continuous Inverses, 
where by Ryo we mean the restriction of the resolvent that A’ is bounded relative to A®, and that the set 
operator for A°® to the part of our Hilbert space p is complete in AH ’ For convenience we will 
orthogonal to 4 The eqs 20) become assume that A®° is positive definite: the modifications 
for A bounded below are straightforward 


" The yur ose of the yroof is to show that the inverse 
S a p R, A’ p r A’p t S a;\i Ap (), . p | — St = ar } 
cud é i f f operators converge uniformly to A From this 

the convergence of the eigenvalues follows immedi- 
® kr) 26) ately [5, p. 372]. 
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Since the inverse operators (A*)~' form a decreas- 
ing sequence of symmetric operators bounded below 
by A@', the sequence converges strongly to some 
symmetric operator [5, p. 263]. The strong con- 
vergence to A™' is guaranteed by the inequalities 
that hold for any ¢t in % 


{(A*) 1_ 4 "Nt ‘ 
<|\(A°)~ 


< —A~*)3)\t)|2([((A*)-'—A-']t,t), (30) 


and 


({((A*)-'—A™'It,t)?<[(UI—P) A“, —P*)A“'t] 
-(A’(A*)—'t,(A*)—'t). (31). 


The first of these inequalities follows by using the 
Schwarz inequality in the form generated by (A‘)~' 
A~' and then using the fact that (A*)~'<(A®)"'. 
The second is obtained by writing (A*)~'—A™' 
(A*)'A’U—P*)A', and then using the Schwarz 
inequality in the form generated by A’. The first 
factor on the right in (31) vanishes with increasing 
k since the set {p,} is complete in #’. The second 
remains bounded since A’ is bounded relative to A’. 
The strong convergence of (A*)~' to A is thus in- 
sured. However, since each operator (A*)~' as well 
as A~' is completely continuous and symmetric, the 
monotonic strong convergence implies uniform con- 
vergency by an easily seuved analog of Dini’s 
theorem. 

Remark: We observe that if. the span of the 
vectors {p;}f.; contains some eigenvectors {u,, U, 

» u,,} of A, (G<k), then A’ 

these u’s so that d,,%,, . . . , A, are eigenvalues of 
A* as well as A, and they have the corresponding 
eigenvectors U,,, U,,, .. . , Uy», respectively. In- 
vestigations based on this property lead to state- 
ments of the following type: 

THeoreM: Let A be an operator with a purely 
point spectrum and let \, be an eigenvalue of A that 
is less than A, then there exists an intermediate 


problem (of sufficiently high order) which has 
M= kX, and uk—w«, for nw—1,2, . .. , ». 


3. Truncation of the Base Operator 


The principal difficulties in the computation of 
approximating eigenvalues and eigenvectors by the 
method just described arise from the fact that the 
resolvent operator fy for the base problem rarely 
is known in closed form. In most cases it can be 
expressed only by use of infinite sums or integrals; 
i.e., in our cases by 


R \ (w,u?)ue "°° ) =6AE®w 

1\W : 

oat Ny r Jin~— Oo d 
We overcome this difficulty by changing the base 
problem so that it has a simple resolvent and still 
provides lower bounds by the method of inter- 
mediate problems. These bounds are computed 
from equations identical with those given earlier ex- 


A on the span of | 


cept that the resolvent is expressed by a finite sum 
with a remainder in closed form. 

Let us assume that A® has infinitely many eigen- 
values less than )e (the modifications for the cases 
in which there are only a finite number will be im 
mediate). Define the truncated operators A by the 
equations 


At °— A° RY +Mii(I1—E! 


IN \ “a 


and denote the corresponding resolutions of the 
identity by kK. [tis clear that each operator A 
is bounded and symmetric, has the eigenvalues and 
eigenvectors \°, vu? for vy=1,2, . . ., 7, and also has 
MZ, as an eigenvalue of infinite multiplicity. The 
resolutions of the identity may be written down im- 
mediately. They are 


An operator A‘’ ts called a truncation of A® of order 7 
These operators satisfy the fundamental inequalities 


ACCA vg (a=1,2,...), 34) 


since from the spectral theorem it follows that 


({ Att! °— Ahlan AQ o—A d( Fyu,u)>0 


for each win % and that 


({A°—A**|u,u)= (A—XA 


for each win Yo. 
The new intermediate operators A 
A*°® are defined by 


basec on 


AtkF— Ato A’pr G,k=1,2, ...). 36) 


These are bounded, symmetric, and are monotoni 
cally increasing in 7 and k; i.e 


OS) 


The montonicity in the first index follows from (34) 
and in the second index from (13). The inequalities 
(37) and (38) imply the parallel inequalities for the 
corresponding eigenvalues, 


p=—1,2,..-), (S59) 
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where by A) * we mean the vth eigenvalue of A A'* corresponding to °, then uw is given by the 
Since A and A‘® differ by a completely con equation 

tinuous operator, the limit point X! persists and 

is the only limit point of the spectrum of the inter " =. 2 1{’n S° » (44) 

mediate operator A as the order k of the operator = | ye in E =k+ 

is increased [5, p 367] Further, the inequalities 

39) show that for any & there are at most 7 eigen- | why re the constants a; must satisfy the equations 


values of A‘* less than A 
The spectrum of each operator A'* can be deter 


mined easily by the usual procedure for intermediate hen Hi PT hi 1 pP 1 p 
problems. The eigenvalue equation for <A‘* is 
stated by + » att Ap {). } 4 kf) (45) 
1'*u—du—0 
ind 
or equivalently (in the notation of the preceding 
section dja A’p / ), (7 k- 1 k+2 Th l). (46) 
a | , ri . 
A‘ "u—hu — Ss ja A P); The linear independence of the eigenvectors arising 


from linearly independent sets of a’s mia be demon- 
strated as in the preceding section 
SS a,(p,;,A’p u,A’p,), (i=1,2,...,4 The multiplicity of A2 as an eigenvalue of A‘* is 
thus just the nullity of the matrix 


If w is an eigenvector corresponding to a value A 


that is not in the spectrum of A‘t®, then « may be Py+RyoA’p,,A’p,) ,A’p 
written (47) 
lp ft 0) 
u — Sa R\A’p,;, (40 

We have found all of the eigenvalues of A‘* that 
ure not equal to X A The complete continuity of 
where RY is the resolvent operator of A We A’'P implies that X°.,. the only limit point of the 
note here that for any vector v in #% Rip is given spectrum of A®°. remains as the onlv limit point of 
bv the closed expression, the spectrum of A Thus the total spectrum of 

A Is determined The spectre al subspace of A 

R wh (0, uy) u | i en, 7 associated with NX, is just the orthogonal comple 
*) —— y \ Tt \ r a y Uy) U | ment of the span of all the eigenvectors already 

PS determined for the values of \ not equal to AY, 
+1 Thus. as stated at the beginning of this section. 
the method of truneation justifies the replacement 
The eqs (22) become of the base problem resolvent by a suitable approx- 
imation in closed form. The necessary computa- 
~~ > " . a -.. : tions can be carried out by routine techniques. 

2 a@(pyt+hy A’p,, A’p 0, (=1,2, A The convergence of the eigenvalues of the trun- 
cated intermediate operator A'* to those of A* can 
(42) | be shown under the assumption that A® and A have 
completely continuous inverses We demonstrate 
As a consequence, the necessary and sufficient condi- | the uniform convergence of (A‘*)~' to (A*)~! for 
tion for such a A to be an eigenvalue is that it be a each fixed j For convenience we will assume that 
root of the equation A° is positive definite. First we note that the 
operators (A ' converge in norm to (A°)7! as 7 
Pyt Ry A’p,, A'p,)|=0 13 becomes infinite. This follows directly from the 


spectral theorem, for we have 


The multiplicity of each such eigenvalue is the nullity . 
of the coefficient matrix of eq (42) The eigen- (A dk. 
vectors of A‘* are determined from the solutions ‘ 
of eq }) by eq 10) 
Now suppose that 2 is an eigenvalue of A that anc 
is less than A and that the characteristic sub- 
space 4 of A corresponding to A’ is spanned by A h~'dE°4 dE®. 
i If uw is to be an eigenvector of e A Jd 
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so that 


(A%)-1—(A*®) =| (j-n-)4 zo. 


Consequently we have 


, l 
|| (A%)-'—(A*)-" || <S — (48) 
: —~ NG4t 
Since \f,, becomes infinite with 7, the norm = con- 
vergence of (A‘°)~' is clear. This implies that 


(A**)~' converges in norm to (A*)~'. The proof 
goes as follows. We start from the simple identity 
At*®— A*= A‘°— A? 


From this it follows immediately that 





At.*|(A*)-'—(A**)-"|A®= A" (A?) (A**)-"|A° 
and thus 
(A*)-'—(A*%*)-!'=(At*)-' 44° 

[(A°)-!— (A*°)—-"|49( AS) (49) 


By writing (A**)"'A*® as I[—(A**)"'A’P* it is evi- 

dent that (A‘**)~'A*® is an operator that can be 

F aued Mal be 4 

=X 

Since A®(A*)~' is also bounded for each fixed k, the 

convergence in norm of (A‘*)~' to (A*)~' follows 
directly from (48) and (49). 

The uniform convergence of (A**)~' to (A*)7! 
and that of (A*)~' to A~' demonstrated in section 2 
are sufficient to insure that the eigenvalues of A 
converge to those of A as 7 and k become large. 


bounded independently of 7 since 


4. Truncation by Choice of Elements 


In this section we discuss an extension of another 
procedure [2] for truncating the resolvent 2, that 
appears in section 2. This can be accomplished by 


choosing the elements p, in suitable ways instead of 


changing the base problem operator as has been 
done in previous sections. 


We assume that the vectors { p,, po, .. ., py) which 
determine A* have been chosen so that 
N 
A’p, >> Bp u?, (s=1, 2 . k). (50) 
v= 


That is, it is assumed that each of quantities 
A’p(i=1, 2,..., k) can be expressed in terms of a 
finite number N of eigenvectors of A®°. The number 
N will of course depend on k. The quantities 
(R,A’p;,A’p,;) which appear in section 2 will have 
the simple expressions 


. - 
SO By B, 


; (51) 
—S \"—» . 


(Ry. 1 "Diy | ‘p ) 


and thus are effectively truncated by the choice of 
elements (50). The computations for the 
values and eigenvectors may be carried out directly 
from the equations of section 2 

In the present however, it is simpler to 
proceed directly. We have already shown that the 
expression A*u can be written 


eigen- 


cause, 


Atu=A®ut+ >) (u,A’p,)b),A’p,, 


where, as before, the quantities 6,, are the elements 
of the matrix inverse to that having (A’p,,p,;) as 


elements Using (50) this expression becomes 


\ 


Atu=Aut > 


We can now obtain all of the eigenvalues and eigen- 
by inspection First 
that are orthogonal to every 
If wis such a vector, then (52) shows 
that A*u—dAu— A°u— dru ( leary no such a vector 
be an eigenvector of A*. Further, if uw is an 
'with an eigenvalue X\ 
then by 52) we have 

Hence such 
with the same eigen 


vectors of A’ consider those 
vectors 1D ‘S4 


vector of A° 


even 


can 
eigenvector u° of A 
not used in the sums in (50), 
A*ue— dru? =A°u rou? —0 
vector 1s an eigenvector of A’ 
value as for A 


and Is 


each 


It remains to consider those vectors uo which 
are linear combinations of the vectors u(y l 
2 .. N) used in the sums in (50 There are 


exactly N linearly independent vectors of this kind 


that are eigenvectors of A’ These are determined 


as follows. We write was 


\ 
and insert this in the eigenvalue equation 
A*u—dru=0, where A* is given by 52). The 
linear independence of the vectors wu leads to the 


equiy alent algebraic system, 


Consequently the determination of the intermediate 
eigenvalues and eigenvectors depends only on the 
solution of a matrix eigenvalue problem 

This procedure can be especially useful when the 
multiplicative operator, as for 


operator A’ is a 
the 


example in the spheriodal Wave equation.” In 
prolate case the eigenvalue equation 1s stated by 





which A is essentially self-adjoint operator 


in 
riven by 


on the set of all functions continuous on 1] 
that are twice continuously differentiable on 
l, 1) and for which Au belongs to Y l, | 
Here mis a non-negative mtegel and cia real 


constant \ base operator 1s obtained by omitting 
The solutions of the base 
nSssor inted Legendre The 


and the recursion relation 


the positive teem A’ ( 
problem are functions 


choice p / PP 


m)P 


allow as a linear combination of 


u , w 


me to express 
Ue. Wf ie and / 


this sort can be found easily 


Othe examples of 
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1. Introduction 


M.. Stone obtained the fol- 
theorems of 


In a recent paper 1] 
lowing result in reproving 
Khrenfeld [1]: ¢f F and G are positive definite n-square 
Tlermitian the) 


certain 
matrices 
d( I+ AF) >d(1+-AG l 


Tlern itian matrices A 
semidefinite Hermitian 
henceforth A >0O 
positive semi- 


for all positive definite 7 
j fk G 


square 
if and only 
Here determinant 
>) will signify that A is positive 
definite) Hermitian. /F'>G (Ff >@) will mean F—G 
-O (fF G >0 Note that since A Oif and only if 
A 0 the condition (1) is the same as saving 


S positive 


d denotes and 


d A I > «l A G 2 


for all A>O It is in the form (2) that we investi- 
cute W h il happens when we replace d by some othe 
elementary symmetric function of the characteristic 
To be specific, suppose Kk I ap ] 


roots . ’ 
: elementary svmmetric function of the 


denotes the 7 


indicated variables 


fixed r, 1<r<n, let (A) denote /,(\,(A), 
(A where A,(A), 2 :. _ nm, are the 
characteristic roots of A It the A,CA) are real we 


Kor il 


will choose our notation so that A,(4) : 

(A Then the problem we pose Is to find con- 
ditions on the characteristic roots of /’ and G@ such 
that 


A+G) for all A>0O. 


We hav c 


THEOREM 1. Assume F'>0, G>O0, 


10, 1961 


‘CSSal ind Fic conditions f 1-/ 
nditior ! ms of the compa n of 
eris fA Y and B Ya \ i : 
lf 
| r,, (F d F) | LilG d »(G@)), 
(3) 
p ) } ] 
it follows that 
/ | I I | Gr) for all A 0) 4) 


3) become S 


In case G /, 


which is both nece ssary and suffic ent for 


BE(A+ F)>E(A+] for all A>O (6) 


We remark that for » ? the second 
theorem | simply becomes: \,(/ >1 «af and only a} 
d(A+F)>d(A+T]) for all A>O. Now im (1) 
multiply both sides by dLA~'R*R) where R is a non- 
singular matrix satisfving R*GR=/], R*FR=K 
and we find that 


part of 


d(R*A 'R K > R*A R /). 

Now as A runs over all positive definite matrices so 
R*A R and from the above remark we con- 
clude that A, (C4 > | The characteristic roots of K 
are just the characeristic roots of G-'F and thus 


does 


G-'2FG-"?> I, h>G, 


and the result in [4] follows. A somewhat stronger 
multiplicative analogue of theorem 1 is available 

THEOREM 2. Assume F>0. G>O. and 1<r<n 
Then 


Ie AF) >I (AG) for all A>OoO 7 


113 








r 
if and only if II A,(F'@) <1. (8) 

j=l 
In order to prove these results we use some results 
concerning compound matrices and Grassmann prod- 
ucts [2]. We give the coordinate definition of these 
items and list several of their properties. If sv, 
., Z, are vectors in the unitary space of n-tuples, 
S.==(Ga, ...,Fa), 1Sr<sn, then SA... AZ,is the 


n 
}tuple whose coordinates are the r-square sub- 
, 


determinants of the rn matrix (7,,), i=1,. . ., 7, 
j=l, .. ., n arranged in lexicographic order. If A 


A) is the ( )-square 


matrix whose entries are the r-square subdetermi- 
nants of A arranged in doubly lexicographic order 
according to the row ~ column ie es of A. That 
is, if 1 <i; -<n and 1< ae den 
are two increasing eo of r helm then the (i, 
er * eS jr) element of (CA) is the deter- 
minant of the matrix A[7,, . . ., ils, : j,| whose 
s,t) entry isa,,. Ift=j,, s=1,...,7 we denote 
the corresponding principal submatrix by A[#,, 


is an n-square matrix, then 


?,}. ge ep rege eed Aer A. sd tmAZa 
Agz,;if (, ) is the usual unitary inner product in 
leaanaee of m-tuples then (4, a... AZ, Yr aA... A Yr) 


Silt, ¥)}, 1*,37=1, . mar 

We remark that the condition theorem | 
will not imply that k>G in the case r<—n. For even 
if G=J, we can easily construct F' so that A,(/)<1 
and yet (3) holds for p=0, ..., r—1. However, 
there is a direct generalization of Stone’s result to 
the compound matrix. 

Tuoeorem 3. Jf F>0, G 


(3) of 


0 and then 


l<r<n, 
1G) for all A>O 
| and only if F>G. (9) 


Once again, r=n is precisely the result in [4]. 


2. Proofs 


We prove theorem 1. Let A=UXU*, U unitary, 
X=diag (x, . . r,). Then £,(A+F)—E,(A+ G) 
E,(U*AU+U *FU) E(U*AU+ U*GU)=EVX 
kK)—E,(X+ HA) aise re kK=U*FU>0, H=U*GU 
>0. Let 
g(r, . rn. UN=E,(X+-K)—E,(X4+7)). 
Then 
@ (21, .. +) 22: U) ae : > — d(diag(r.,,.. -,%a,) 
+K [au,...,a,]|—d|[ diag (z.,, . . ., La.) 


+H [ay, .. ., a,]] (10) 
We see from (10) that ¢ is a polynomial in z;, 

x, of degree r in which the nonconstant terms form 
a multilinear polynomial. Let 


0’ 


Oa, ... O2;, t,:U) 


- 7 


Pip reve ty (Ly + + 09 Myr UV) = 


and 


Then 


where the term corresponding to p=0 Is just g(0 
,e € Let w denote the 


i i,. Unless wis a subset of a Ge a, 


Increasing sequence 


the partial derivative satisfies 


O/ ; 
d diag (27, ,* » La K |a ce 
O@, ...@a 
)? 
d\idiag a » we Tlie ‘ , 4 () 
If wis a subset of a a let a a Lye 
the ordered complementary set of w in a a 


In this case the derivative in (12) has the value 


d{diagw (2.0, - K [a¥’, , a 

didiag Ses’. ae /] a. a , > 
Setting 7 Bern J 0 in (13) we have 
Q p hy d K Gn. 4 

. l<a@ 7,“ 
dH «ct . ° a> : * 14) 

where » Dey indicates that the summation is taken only 
over those a na a, containing w ‘ l us 
asubset. Returning to (11) we set a ! fand 
all other a 0. Then (0... .. ¢. - oe S| 
to + L,-1(t) where L,_; isa polynomial of degree 


at most p—l in ¢ since (7, , a U)>0 we 


conclude, by letting ¢t increase, that 6°... >0. Con 
versely, if ¢!,...,.>0 and (0, . , 0, U)>0, then 
OZ... « oo Zar ) for all non-negative a;, J 

Let J , Is hen the set complement: iry to 
7) igs 2 in l, ——eE © Let t, 3 l, . . be 


the unit n-tuple with 1 in position j, elsewhere 


We may rewrite (14) as 
o} i p dy dk te a eae dH 01. oO 
- log ( # (Kes a.. A é, 5» &«,A Ag Pe 
—(C,_,(H)ee,.A... AG 9 fa, Ae AG 
bby 44 p(P )ug A... A Ue yy le,A- +A Ue) 
—S>’'(C, p(G@)Us A... A Us 
Us A Ale ) (15) 





where » tn indicates that the summation extends over This last matrix has the same roots as the Hermitian 
precisely those Increasing sequence @ J matrix ¢ A I ( A G and hence (9 is 
which are subsets of J Also r lie equivalent to 

] | mois an orthonormal set of vectors (recall 

that @° is unitary We use an extremal result in C'(A+H4 C'(A+G for all A>0O 

}: theorem 1, p. 525] to conclude from (15) that Now 


/ N(G \ G () CG AG LCG CG 1G 
hn case G / we have from 15) again that G KG ‘ G AG f) 
d > ( F yu, a A Thus (9) is equivalent to 


We A Au ('(A+M)—CAA+1)>0 for all A>O 18) 


for all choices of sets of n p orthonormal vectors where /J/=G 2HIG 

u l , and another application of the above 

cited extremal result completes the proof By a unitary congruence we may assume //=diag 
To proceed to the proof of theorem 2, choose a h h,) and by setting A= diag (a7,, ,2z,)>0 


nonsingular R such that F=RR*, G=RDR* we see that 

DD diag A / (7 \ / G Let P be an 

arbitrary nonsingular matrix and since any A >0O ts r 

of the form (P?R gl 3 we have that (7) 1s II (4 hi, at ' 


equivalent to 


/ PR-")* (PR) RR*)> I PR-')* (PR-') RDR* for any non-negative numbers 4 J This 


clearly implies that each A ] t a , = 


(PP i P*PD Thus, 0>H—I=G-“3FG@ [—=F-—@G 


IA) > FAD) for all A>0O 16 Conversely suppose F—G>0. Then //>TJ and if 
we set B=A+J]>0 we would like to conclude that 

In (16 replace A by VNV™*. V unitarv, V=—diag 
J 7 () to obtain f A // ( RB T] / >/ B ( A ] 19) 
K(X Ie NH H7T—VFDv But (19 


is equivalent to 
Oo! 


ire X) U— CD) =0 Mi C+K)>C(D), K=B-"*(H—DB-"2, (20) 
where / is the square identity matrix. It is not 
diffieult to check that (17) holds for all non-negative 
diagonal 4 if and only il every diagonal element ol 


/ ( i] Is non-necatlive That is ( / diag k,. - : I >f /). 


After a unitary congruence 20) reduces to 


| ( V*DV )ea.A A@e. fa.A A¢ >) where k, >) are the characteristic roots of K The 
prool is complete 
must hold for all V and all l<a aS? But 
this is precisely equivalent to 
Citta A. he. a... dn. <1 2. References 
for all orthonormal u,.. , Ge As in the proof 1] 8. Ehrenfeld, Complete class theorems in experimental 
: oA . design, Proc. of the Third Berkeley Symposium on 
of theorem 1 we have finally that é holds if and Mathematical Statistics and Probability 1, 57 (1954 
only if (8S) does 1955), Univ. of Calif. Press, Berkeley, 1956; Math 
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lo prove theorem 0 it will be convenient to let 21 W. Graeub, Lineare Algebra. cl 5 (Springer-Verlag 
denote the relation of Hermitian congruence Chen Berlin, 1958 ~ 
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Some Properties of the Empirical Distribution 
Function of a Random Process 
M. M. Siddiqui 


January 24, 1961 


Let \ i i continuous-time continuous-in-the-mean real strictly stationary 
random process If c is a given number, denote by pr(z) the proportion of the time A 
exceeds O<1(<T Phe covariance of pr(z and pris s Obtained An approximate 
solution for the variance of i sample quantile Is ALSO give! Ihe general results are specia 

ed for Gaussiat ind x processes 


1. Introduction 


The proportion-of-time distribution plays a very important role in the stitistical analysis 


of time series data. The practice is to plot p(s), the proportion of the time the record, X(t 


exceeds r, against J 6... 1,2] In veneral if P { t Is a wide-sense stationary process and a 
sample over O0<t < T is available, an estimate, Pris of the probability that YV(¢ > ris obtained 
from the sample. In this paper we will derive the covariance of pr(z,) and pr(s An ap- 


proximate solution for the variance of a sample quantile will also be given. The general results 
will be specialized for Gaussian and x” processes The Ravleigh process Is a special case of the 
x’ process; however, due to its wide applications, the result for this particular case will also be 


stated 


2. Some Properties of Time Averages 


Let NX, (t X; t)}, teR?, be a wide-sense Stationary, real, continuous parameter, con 


tinuous-in-tie-mean, vector process t will be called time and R will be taken to be a real 


line 
For i=1, 2. and for all ¢ and s, let 
BX ,(t)=m 
EX (t)X,(t+s)—n Ri(s), 2 
EX, (t).Xo(t+s8)—mmo=Ry»(s 
where /.XY denotes the mathematical expectation of the random variable Y. It will be assumed 


that all quantities in (2.1) are finite 

When statements are made which apply to either of the processes, O1 when the properties 
of a single process are considered, the process will be referred to as X t ; and the subscripts 
will be omitted, e.g., A.N(t m BN (t)N(t+s n Ris 

Let T be a positive number, \ a positive integer, d= 7/N, and denote X,(kd) by X,(k 
Ritkd) by Rh ete., since there will be little danger of confusion X,(k), Notk k—0O. 

a 2 _is, then, a wide-sense stationary, real, discrete-parameter, vector process 

We may consider the discrete-time process and extend the results to the continuous-time 
process by taking the limit as V->, 7 fixed; or, we may consider the continuous-time process 
and obtain results for the discrete-time process by replacing integrals by appropriate summa- 
tions. In general, we will pursue the latter procedure 








Define 


Let t,, t, be arbitrary numbers, and 7), 7; arbitrary positive numbers 


It will be assumed that the subscripts on Y are so chosen that 7,2 7). 
For the discrete-time case, we take k,, k, arbitrary integers; N,, N, arbitrary positive in- 


tegers, choosing the subscripts on Y so that V,2.N,; and define 


= 1 Nut 
Y iw, (k;) V pa =1,2 Z.2) 
Now 
EX,, (7;) Mn ee. 2.3 
and 
cov{Xir,(t:), Xor,(t,)} =(Ti M2) EF ~ LXi(t’)— mi] [X2(t) — m,| dt dt’ 
*7 *T, » 
(T,T.) R(t—t’ dt dt’ 2.4) 
Jt vt 
The interchange of operations of expectation and integration in (2.3) and (2.4) is justified 


as HN? (t o, 4=1,2. 
We now make the transformation u=t—t’, 
unity. The limits of u and v are given by 


t’ The Jacobian of the transformation is 


t.—t:—T,s ust.—t,+T, 
max (¢;,t.—u)Svs min (7,+1¢,, 7.+t.—w). 


Hence (2.4) becomes 
*T)+8 
(7, T;)= (min (7+, 7.+t,—u)—max (¢), to—w)) Ry (u)du, 


e T; T 


where s=t,—t,. Now 
7, +1t,, when uSs+T,—T 


min (7,+¢,, 7.+t.—u) 
T.+t.—u, otherwise; 


t,.—u, when uSs, 5 
max (t),f.—wU) ; (2.5) 
t,, otherwise. 


Replacing ¢, by t, we obtain, for arbitrary ¢ and s, 


rs 


= oo s °T,-1 
cov {X,7,(t), Xor,(t+s)}=(1T2) i (T,—s+ wu R,.(u)du-+ T; Ri.(ujdu 
; . . T,+8 8s 


T 4. 
+ (T: | o— ay) f? 9 | udu | (2.6) 
7 Tl, +8 


Remark 1. Setting X,(t)=X,(t)= X(¢) in (2.6), we obtain 


—_— = "s 
COV {Xr,(t), X7, (t +-s)}=(7;T)) ‘| (T,—s+u)R(u)du 
—T,;+8 


_ (T,-T,+8 *T,+8 - 7 
+1] Rowdut | (T, +-s—u)R(u)du } (2.7) 
s r,—T\+8 4 
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Remark 2. If 7T,=7 / 
cov {X7(t), Xr(t+s ] | T—s+u)R(u)du4 ¥ T+s—wR(udul (2.8 
Remark 3 Setting s Qin (2.8 we have 


var Xe(Q=T'f (17) Rwdu 2.9 


From (2.3) to (2.9) we deduce that, for fixed 7) and 7 A t A l is a wide-sense 
stationary continuous-time and contimuous-in-the-mean process 
For the discrete-time case, we will stute only the results corresponding to (2.6 and 2.9 


When translating a sum of integrals, such as 


into a Riemann sum, we should be careful to note that, whereas 6 appears twice in the limits 
ol the inte erals tiie corresponding quantity will appeal only once in the sum ol two summations 


Thus, for arbitrary integers k and 


cov {Xv (bE), Xev (k-+s VA > \ uw) Rou 


/ Z 
Similarly 
— RO 2 X-! l ; 
var Ay(/ { rae ) Ru 2.9 
\ y \ 
The results (2.9) and (2.9’) are well known [3, p. 80]. 


3. Empirical Distribution Function 


It will be assumed further that {.V(¢)) is a strictly stationary process so that the probability 
distribution of .YV(¢ h Vi / is the same as that of YV(? Y(f lol every 
selection of f f and / 


Define 


. 
I f-rjydt 5.2 
| | 
is the proportion of the time s Is exceeded by Vi?) in the time interval (0. 7 i. ] r) is the 
La hy syue measure of the set S / Q </-< ] We note the followin Vv 
j | for all ¢; hence p l 
2 pel UO tor a | f- hence p ) 


3 lf ry, then X(t) 2a, implies V(t) 24 


Ss is, theretore, a subset ols and hence 


Properties 1) to ‘ imply that g-(a | pr(r) is a distribution function dr\4 will be called 
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the empirical univariate distribution function of the process 


N(t We will, however, consider 
pr(#) rather than qr(z). 
Suppressing « from p(t; «) for the time being, we have 
Pr[p(tith)=1,..., p(et+h)=1, pit h) pt h 0 

Pr(X(ti+h)227,..., N(t+h)27, X(t hy<a N(t,+h)<a 

PeiAt,)zeez...., Alb)ez, X(t ! NX(t } 

Prip(t,)=1,..., plts)=1, pltes:)=0, pit 0) 3.3 
for every selection of k,n, t,, .. . . tt, ...,t,, and A. Hence, for a given «, (p(t: 7)) isa 
strictly stationary binomial random process 


values, Zi, _ = pit, 2), 


In fact, it can be shown that, for 





a given set ol 
‘ . pit; I, Is a strictly stationary vector process 
Writing 
P(x) =Pr(X(t)= 2), 3.4) 
we have 
Ep’ (t; x) P(z), 1 1,2 ) 
var p(t; r)=P(2)[1—P(2)), 
> (3.5) 
A(s: r)=cov p(t; r), p t-+-s: zr) 
Pri X(0)>27, X(s)22 P2(2 J 
Let 2 =27,; then 
A(s: 2 ro) =cov p(t; Ty , pt 8; Do) 
Ep(t; x;)p(t+s; x, P(2,)P(22) 
Pr{ X(t)2>”, X(t+s) +, P(2,)P (2, 
Pr{ X(0)=>2,, X(s) Pe P(2,)P(a). (3.0) 
We note the following: 
E’pr(r)=P(2) 
A(s; z, r)=A(s; 2), 
A(0; 21, %)=P(z.)(1—P(2 
Using (2.6) with s=0, 7,;=—T,=—T, we have 
+ 


1 ¢? ‘ 
O21, Zo) CoV Ie D1), Drl Le) = |— — A(s: a”. 2 ds, 
I 2 Pr Pr\t2 7 a 7 ) 


C(2, z)=var pr(z) 


For the discrete-time case, 





. N 5 
Pr(2)=N-! SS p(k; 2), 
k=1 
: \ 
COV { Py(%1), Pw(22) } = N~'P(2,)[1—P(a,)|+2N- S53 —s/N)A(s; 14, 22), > (3.7’) 
s=l 
N l 
var py(r)=N"'P(2)[1—P(z) 2N-! 3) (1—8/N)A(s; 2) 

¢= | 4 








To estimate A(s: 4 j and A Rid define 


lif VC? + V(t . +” 
p - \ J I 
0 otherwise 


aha use 


Lo estimate P1 X f sd, X(t \ = 


4. Variance of a Quantile 


Griven J p j Is a random variable In the preceding section we obtained the variance of 
_ and the covariance of Pr\a and py J In this section we will obtain an approximation 
to the variance of a quantile of the empiri al distribution function We suppose that a number 
q ()} q l. is sper ified If [ denotes the number suc h that 
p { , +. 
then 0” is a random variable. Let & be the number such that 
Te , 9») 
P(t) =, 1.2 
1. ¢ fis the qth quantile of the univariate distribution function of the process LV (t)] 
In what follows the primes will denote differentiation with respect to the indicated argu- 
ment Furthermore, it will be assumed that P? é and Pris have derivatives. 
Now 
pP be pP U] pP Ps [ jp = 


where= means ‘approximately equal to:”’ hence, 


Now 
Ep',(¢)=E lim ! 


theretore 
Coy [ P17 bs = pr t kad 
As a first approximation we suppose that 
sO that 
Coy Up t ) 1.5 


As a further approximation suppose that U7 and p>(é) are independent. Writing dy=y—E£y, 
where y is a random variable, we have 
bprl&) = — prli él 


var prt =varl Ep; (i 4.6) 
Now 
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(ext€t®) Bett)) (ere t®) Pes?) 


h,k-0 \ h 


Ep? (¢)=E lim 
E lim (hk)~"[6pr(E+h) —bprl(E) +P (E+h) —P(E) | lbp r(E 
hk, k-0 


P’(&)+ lim (hk) "(CE +h, E+) '(¢ 
h,k 


~~) 


Of:0f § 
— : Ort (1, é ™ 
var U7 =| P’*(é)+ ee" ei =',=¢ var prlé), 4.7) 
OF OE J 
and var pr(é) is given by (3.7). The result (4.7) assumes the existence of various derivatives 


In applying this result it is necessary to look into various approxi 
1.7 The approximation (4.5) seems fairly 


and pr(é) needs careful examination 


appearing in the equation. 
mating assumptions which are used in arriving at 
reasonable; but the assumption of independence of [ 
in each individual case. 


From (3.7 


02" (2,. x 1 >? *7 P | 1 sl\ 02A(s 
se Se a | (1-F) AG, n)ds= 7) (1-7) Ls 
Or, O07, 1 O72, O27, e T 1 py I / Or, O07 
in case the interchange of the integration and the partial differentiation is justified If 
, OA 
ds 
7 OX, O22 
then 
Ff | ] 
cS) 0 | ap)? 
OL, OL, ] 
and therefore 
var U’ =var p7(é)/P(E)|1+0(T 1.9 
5. Applications 
5.1. Gaussian Process 
If {.XV(t)} is Gaussian, write 
R(0)=—e V(t X(t ” t Rit 
It is known [4, pp 355-6! that 
Pr{¥(Q2m, V(t+s)2y)=(2e p | exp { > } 7 
. H uy )HT f 
Ply P(y +>) p aa 
Ply) P(y2)+ dS 07(s)t (yr ly 5.1 
Here 
P(y) 2n é d 





H I is the Hermite polynomial ol degree r given by 
H, (x)= d/dr)' « ! y-t4— a rd 


and 


is known as the tetrachoric function of ordet Ps and its tables are available oO}. 


Thus 
a I m J m ) 
Aesa,2)—=SseKor (22) 22") ] 
A(s; 1) =) p'(s)73(“—— ) | 
J 
nserting these values in (3.7), we have, since p s p(s 


Differentiating the first equation in (5.3) with respect to 2, and 2, we get the formal 


expression 


O- f _ : /T) OG 1)p*(s) 73,,(5 m\ ds 5 4) 


If the integral converges, we can evaluate it, and using (4.7) obtain the variance of the sample 


quantile 


A spec al case If Mi, 
Varprl? 7 | ge > s) rr? (O)d 
Now 
) 0, oe 
| 7)! 
T (0) — , 7=0, 1 

) ) ? . 

a7 =#\J «] 1)! 
Hence 

l »7)! 

Var) - ——* 1 p s) ds 5.9) 
rl 2 AST i) 








If L’ is the sample median, 


arc 
Or; Or» 71 


in case the integral converges. 


Hence the series (5.5) will alwa 
integral in (5.6) is that 


exists for some 6,0<—)<. 1 


vs converge. 


then =m, and 


Zz» 


7 (rot) { (1 s/T’) 


°7 
(xo? T) (1—s/T) 


We note that, for large value of / 


(27)! | 
| — 
»)e 


(7!)? \7)) 


A sufficient condition for the convergence of the 


prs 


lim s?{1 


s—>0 


limit is nonzero for some } 





ie >1 then the integral diverges. 
Exam ple 1. If p(s)\=p(—s)=e-“,A>0,s=0, 
7 ; 
(1 s Te Mis — 5. yp tO T é / 1,2 
: Le ae | 

) Therefore 
/ l (27)! l a (27)! 

var pr(m) = >. mm La _ =,” LO(T- 4 
) adi AT = p21( 9!)2(27 +1 wr? 7? 45 2?/( 9!)?2(274+-1 
. = ().35/(AT) —0.33/(AT)?+ O( Te 

Also 
2 >A\I 9.\9 
o* Soa (<9)! i> 47 Lt” T 1, 
OxXO2>|, m Awo?T’ Fy 2(7!)7(27) «= -wd*o? T? Fy 2”(9!)7(27 


ln2 


rho?T rho? 


Hence, if U 


is the sample medi 


/n2 


| me 


Exam ple If p\s) a> 


var pPr( mM) 


The integral (5.6) does not cor 


median by this method. 


Let 
proceed somewhat differently. 


Exam ple Ss. p(s) cos 


Pr 


In fac 


(4.7) shows that, in this case, Ep?(£) does not exist so that py 


Y(t)>0. YV(t- 


0.17 - 
O ly 
Per 1 


an, 


rn l 0.17 
ar a4 T = oaT 
7 2ro® rho? Tl = h?0?T? 


2ro* Var pr(m) I 


Var prim 
Q2in2 0.3477! 
AT wy\7T? 


» = similar calculation shows that 


ols é =) 


fail in evaluating the variance of the sample 


0.47a/T—0.35a2/T? 


iverge and we 


t, an examination of the assumptions leading to the equator 


is not differentiable 


¢ 


(7s kh) Instead of using 5.5 to evaluate Var prin wie 
Note that [6, p. 290] 
= 
s)2U = - Cos p\s OS cos p(s)Sr 
2 2 
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Hence 


Now 


{ T » 
j 27h if 29k: 2 L)A e==(), | 
cos‘; 
“(s—2jk), if (2j/—1)h 2ik, j=1,2 


where 


Val ; ] 
l 4 , 
h 3T?+- if 7 is odd. 
i ‘ ) 
Obviously var pr(7 0 whenever 7 0, / even, i.e., whenever T'is a multiple of 2k, the period 
ol the process and also vat p n () as T 


5.3. x*® Process 


be indepe ndent Stationary Craussian processes with the same 


Let  A,(t 1, 2, 


variance, o?, and the same autocorrelation function 


Define 


will be called a x? process as the univariate distribution of Y(t) isa x’ distribution with n 


Vit 


rrees of freedom which has the density 


deg function 


4 220 (nH /2)| "6 lol 0 
0 otherwise 


)(t)} is strictly stationary. 
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It is easy to verily that 








The characteristic function, (u,v), of [Y(t), Y(t+-s ] is easily evaluated to be 


o(u, v= Heir tieveess 


(1 Piu Pip fan | p*(s)) 











° . . , 1p° S)u 
(1—22u) (l1—2av)-”?] 1-4 
1 —2in)(1—2) 
) ; 
| 1)?P'(n/2+-7) 2p(s ui = =) 
= “wR : — : rn ? d.4 
| j=0 T(n/2)P(g+1) (l—2eu)"?*7(1—2iv) "22 
as 
fur fur 
(1—2iu)(1—2iv)| [((1+4u2) (+42 
the binomial expansion is uniformly convergent 
For a full discussion of the distribution of [¥(¢), Y(¢+-s)] the reader is referred to [7]. A 
technique of obtaining the probability density function of [ (4), ¥(t+-s)] will be presented here, 
which is also applicable to other distributions which can be expanded in a series of orthogonal 
polynomials. 
Note the following: 
l é , 
(1) g(z; n/2) ; du, 
: 2. l 7) 
d? l 7 Lu )-é : 
g(z; n/2) du, 
d ar. l iu) 
(2) It is easy to verify from the definition of Laguerre polynomials [8, Ch. 5] that, for a>0, 
jSa l, 
@ . PGj+1)P(a—j 
Qi 2; a) Qz:a ) ~ . 2 
dz: . / 2’T'(a) 
where 
. R(OTI\= 
L\™ (2) > > ( | ( ) ,? 
j—k/ ke! 
is the generalized Laguerre polynomial of degree /. 
Hence, if n22—n=1 may be reduced to the Gaussian case discussed in 5.1—-, by the 
inversion of (5.7) the probability density function, g(2,,2; n/2), of [Y(t), Y(t+-s)] is given by 
[2o(s)|*T'(n/2+-7) d 
. [2p 2 
g(21, 22; n/2) S - a, / g(21;n/2+7 7 n/2+-7) 
—& in/2rgt+l) d D Ga! 
a 
C(7+1)P(n/2 ' 
g(21; n/2)g(z2; n/2) >} J p24(s)] 2)L 2 ».§) 
2 ee j=0 I'(n/2+4-7) 
Now, using the first equation in (5.8), we obtain 
Pr{¥(t)2yn,V(t+s)z2y.}= 9( 21; 22;n/2)dz.d 
e Yi e Yo 
2o(s) . C'(n 24 ) d d 
Py) PYy2)d+ > WyYi;n/2+7) J(Yojn/2+-7) 


1 Cin 2)TQ 1) dy 


P(y:) Pye) + 9(yiin/2+ 1) g(yejn/2+1) SO = 


Lj2!? (y,/2) D3"? (y 2). : (5 9) 
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where 


P(y ( n/2 5.10 

‘rom . ithe 5.4 thie covariance OL p {sy and put as) 
, 2nl »)T(n/2 0 
A(s:y,y jnj2+1)g(yo;n/2+1) > Li"! (y/2) Li? (yo/2 
: F = | 1 / Z 
Inserting this value in (3.7), we obtain an expression for the covariance of pr(y) an 
p / The Varinunce is givel by setting Y y 

A pee al case j 2 gives a X® process with two degrees of freedom, i.e., the Ravleigh 


process We find 


COV | Prlti), PrlY IY 2)I(Yo5 2 >. L y,/2)1 y 2) | (1 mn) Is ) 


The author thanks Dr. E. L. Crow for his helpful discussions during the preparation of 


this papel 
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Another Extension of Heinz’s Inequality 


Marvin Marcus 


December 6, 1960 


A recent result of Heinz gives bounds or 


he two Hermitiar 


associated th G9 bv use of the Grassman! 


Let A 


n-square matrices and let (/ be an arbitrary 7 square 


nl al B be non-negative Hermitian nl hy 


matrix. A recent inequality due to E. Heinz [3] 
states that 7f AA V*QV>0 and B WU*>0, then 
(du Au B'-°¢||, O<p<l 


Here D4 
Hermitian 


is found in 


wranyua 
The most 
1] along with sev eral reterences to 


Ss non-negative recent prool 
ol (1 
other proots and extensions 

veneralization ol | 


In this paper the following 


will be presented 


THEOREM If A *Q>0 ai dB t= >0 and 
u f are any 2k vectors, 1 i } 
the) 
det (Qu det | (A?u,, A?u 
det 5 3 
I A’u B 2 
where O< pl 


follows we use certain elementary 
kth 
which are found in [5 

LL.eEMMA If H>0 and K>0 and Il—K>0 
CH C.K) >0 

PRroot We mav assume A is nonsingular 
hence /7 will be) by the standard continuity 
ment Now /7— K>0ifandonlvif // H 

I—H Kil >. This latter inequality holds i 
andonlv if every eigenvalueof // KH 
Now the eigenvalues of ¢ H AK) are the 
taken & at a time of the eigenvalues of /7- 1K 


In what prop 


erties of the matrix of 4 C.(A 


COLLpO ll | 
they 


and 
argu 


AK)H 


isatmost | 
produc ts 


\lore 


ove! Hl K has non-negative eigenvalues ind hence 

everv eigenvalue of €,(/7-'AK) is bounded by 1 

Now ( HK (4H (,(A) and this latter ma 

trix is similar to | H ( Ky )\¢ /1)) Thus 
C1 (<A) '7CO (AYO >() 

aie 


(iH C.AK)>0 


Kor completeness we next ine lude nvery brief and 
which relies on the fact [4] 


is monotone ol every ordet 


elementary prool ol | 


that oA \ 0< p<, 
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parts of Q This 


>) means simply that 


the bilinear form associated with a matr (J 
is extended to certa 


algebra 


\ scalar 


order n on a \< h il 


for non-negative A ¢ is said to 


whenever // 


function 
be monotone of 


K >) it follows that g HI ¢ K >). where I] 
and K are n-square Hermitian with eigenvalues in 
the interval a<A<b To see l let ( UH be 


the polar lactorization ol @ Hi>o, tl unitary 
Then the hypotheses are equivalent to A IP>o0. 
Bb UH" >) Setting w= 


we compute that 


Qu Hu, w HPD Py 
Hu, H-?w) |< (Hu, u)(H2°-?w, w) 


HT?u, Hu) 


To proceed to the prool ol 2 let w A A U, 
denote the Grassmann (outer product 1 of the 
vectors wu wu Then. bv the lemma, 
O< ¢ | ( *0) ( | ( 6 i) 
‘ ‘ t 
nie 


Hence, applving (1) to C,(@ C'.(A 
in place ol | A, 


and (x B) 


‘ 
and B respec tively, we have 


det (/u;, ( (J)UL A A Uy, A Al 
( | Wi A At 
( B A A 
( 1’), ~ Au 
( B A A 
1?u a a A? u 








This last inequality is an application of the Hadamard 


determinant inequality to (A?u;, A?u,;) and 
{ (B'-?»,, B’-?v,)} and completes the proof. 

If ist< ... <teSnm and 14),< ... <pecn 
are sequences of integers, then Ali, .. . , Bal das - - + Sal 


will denote the k-square submatrix of A, (a;,,), 
st Si aare 2 

Coroutitary 1. Jf A?—Q*Q>0 and B—QQ*>0 
and 0< p<, then 


met Gis,, ... » Jai, ..~ ,tallP< 
eerie, « + « » Balt, » s+ > Bal 
det B2°-? La, ei tah its yi diy so eg J ] 
Fmoor. Let @%=e., %=¢;, *=1, ... f 


where ¢, is the unit vector with 1 in position t, 0 
elsewhere. 

Let A and B have eigenvalues a>. . >a, and 
rr oe. ae. respectively. A bound for the left 
hand member of (2) may be given in terms of these 
eigenvalues as follows. 

COROLLARY 2. Under the conditions of the theorem 


rw 


k 
det {(Qu;, v;)} |< Ia? Bi~? (det { (u,, u,)})' 
1 


(det {(v,,v,)})”%. (3) 
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PRrRoor Note that 


det | (A??u,, u C.(A’)m, 


) 


Applying this to 2) we get (5), 
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Mean Motions In Conditionally Periodic Separable 
Systems 


John P. Vinti 


November 10, 1960 


\ search of the literature failed to disclose inv general statement or proot of i theor 
ily current among dvnamical astronomers Ihe present paper give i proot | 
which states that, inv conditionally periodic separable system the mean fre 
quency n, of any separation coordinate g, is equal to »,= 0a,/0J,. Here a, is the energy 
j the Ath action variabl The proof is carried out for nonsingular Staecke 
so that tis appli ible to ! onpolar orbit of at rein I iteiiite when the potential 
ids to separability 
1. Introduction lem [1]? for which the coordinate frequencies are 
all equal to the constant value » v,=v,, and the 
Conditionally period separable systems are com problem of two centers, which has remained a 
monly illustrated in works on advanced dvnamies by | Curiosity to be found in Charlier’s famous book [2], 
the motion of a particle under the joint action of but without application since 1957 Various po- 


‘les This l¢ ntials have been suggested, by sterne 3], by (Ciar- 
happens to be a very special ease. for which each finkel |4], and by the author 5] for the gravitational 


field of an oblate planet, all of which lead to separa- 


harmonic oscillator forces at right ans 


rectangulat coordinate q has a constant lrequency, 


equal to the corresponding ‘ fundamental lrequency bility and to intermediary satellite orbits with 
Oa,/O0/J,, where a, Is the energy and J the cor- dha La Since the solution for these orbits is 
responding action variable greatly facilitated by knowledge of the mean coordi 
In a more general system of this type, each gen- nate frequencies, it now becomes desirable to have 


a formal and veneral prool of the theorem 

It is here convenient to discuss briefly the general 
plan of the paper, without definitions The difficulty 
in proving the theorem arises only when the funda- 


eralized coordinate g ria have a variable frequency, 
but it uppears to be generally believed nmonyg aly 
namical astronomers that the mean Irequency ol q 
must be equal to v,, if the conditionally periodic 


svstem 1s separable Needing to refer to such a mental frequencies Vi, , VM, ure incommensura ble 
theorem in solving a specifi problen 1 have If wy, W, are the angle variables, the plan is first 
searched the literature but have found no explicit to show that for a conditionally periodic svstem 
statement or prool of if The present paper is an there exist infinitely many values of the time, with 
attempt to furnish such a reference, with a proof ho uppel bound, at which the orbit in w-space passes 
sufficiently general to be applicable to all the sepa- arbitrarily closely to points separated from the 
rable problems that may arise in the gravitational initial point w, (0 wy (0 by integer intervals 
theory of the orbit of a satellite of an oblate planet Aw " h I, y This fact follows directly 
Little of the analysis in this paper can claim to be from a theorem of Dirichlet, which is easy to under- 
really new Much of the pertinent material in the stand and to apply 
literature, however, is discursive, relatively unavail- To convert this result to q-space, it is necessary to 


able, and expressed in notations now unfamiliar to know that q, is a single-valued, continuous, periodic 
most mathematical physicists Some of it is inade- function of the w’s lo show this the author re- 
quate, 1 not incorrect, especially in the treatment of stricts considerations to nonsingular Staeckel SVSs- 


the periodicity of the g’s as functions of the angle | tems, proving that for them each g, is a single-valued 


variables Moreover. none ol it seems to have been differentiable function of r, dq pi and each p, of 
carried out in the shortest and most appropriate wa ‘ 
to prove the theorem in question The present the angle variables Ww), , W With a careful use of 
paper attempts to give a concise and correct treat the single-valuedness, the periodic property then 
ment that will serve this purpose follows 

It is easy to see why such a theorem should have Application of these properties of the qs as func- 
escaped formal statement and prool Physicists tions of the w’s shows that, at the values of the time 
have not been concerned with mean frequencies of | mentioned above, the orbit in g-space then passes 
this kind. Dynamical astronomers have been, but | arbitrarily closely to points where each g, would have 
ordinarily for nonseparable svstems Until 1957 gone through exactly m, eveles The proof of the 
their only separable problems were the Kepler prob- theorem then follows 

| i pi ‘ S.A t e. through the Office ~ 

RK t \ Rese 1 pment ¢ nmal Figure ket licate the lite tur el t the end of tl paper 
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2. Staeckel Systems 


If qi, ..-,G,and py,..., p, denote the generalized 
coordinates and momenta of a dynamical system of 
n degrees of freedom, the system is said to be of the 
Staeckel type [2, 6, 7, 8, 9] if the Hamiltonian /7/ is 
given by 


l 2. : 
H 9 2p Ad(a, - - » In) P?t + VK, ~~ In)s 
A. >0, (k=1, .nm) OL) 
and if there exist functions @ (gJWilg), tg=1, 
n, such that 
A, My det (@,;,), (2) 
VSO viele) Ay, (3) 
c=1 


M, being the cofactor of ¢,.,(g¢,) in the determinant 
det (@,). 


Conditions (2) and (3) are necessary and sufficient 
for the separability of a system with such a 
Hamiltonian. 

If we next define the domain QY of the q’s as the 


totality of real values of Gs, « » for which 
pr 29, k=1, . « W, We may then define a non- 
singular Staeckel Svstem as one for which Y Ye) and 
deilqe), 1, K=1, , hn, exist and are single-valued 


and for which det (¢,,) #0 any where in J If we 
put 

P= (¢,,(q;)) +) 
1 


for the Staeckel matrix, then @~! exists and is single- 


valued anywhere in (; in particular, 
(@-'),,.=A, (F=1,...,%) (5) 


all exist anywhere in &. (This restriction thus rules 
out polar orbits from consideration if the right 
ascension ¢ is one of the coordinates, since A, then 
becomes infinite on the polar axis.) 

The momenta p, are then given by 


Pi Wel Ge) +2 DO deil eden: (k=1, ,n), (6) 


i 


where the a’s are separation constants, a, being the 
energy. (For satellite problems, where n=3, 2a 
and 2a are usually denoted by OG and a3.) 

The Hamilton-Jacobi function W is then given | 


"< 


— 
le 


W=>> | ‘pag, 


2 } yf | - 2, ( dy) 2>% 4, we | “dd, (7) 


! 


where the sign is respectively as dq,=0 


| 3. Conditionally Periodic Staeckel Systems 


We call a Staeckel system conditionally period 
if each coordinate is either rotational or librational 
A coordinate g, is rotational if: (7) it is an angl 
i) with p2?=F;,(q,) there exist positive real numbers 
Cy, and ¢, such that c,= fk q 4 0 for all real 
values of g,, and (177) if 


Note that ¢ 


that the periodicity implied by (S) may reduce to 


0 rules out asvVinptotic motions and 


simple constancy The latter holds. e.g., when 
is the right ascension @ of an artificial satellit Spee 
there does not exist an\ potential, depending on thr 
right uSCeHsilon @, which both leads to separability 
and remains finite on the polar axis For sucl 
rotational coordinate g, either p,><« for all gq, or 
P ec? for all q In either case 

q dqy/p 9 


is a single-valued function of g with derivative 
dr,/dq, existing and differing from zero for all values 
of gq Thus gq is a single-valued differentiable 
function of ¢ 

\ coordinate q; is librational if there exist real 
numbers a,, 6 (’,, and ¢ and a real funetion 
G q such that 


p Gi —4,)(b,—qdGi(q 10.1 
with 
( Gq ( 0 a,sy h 10.2 
and 
a,sq (0 h 10 
q 0) being the initial value ol q If we nvull 


define 7, by (9), then 
vA qs) q a,)(b q G q dq, 1] 


where the sign is accordingly us dq —V), respec 


tively Then 


where the unilormizing variable ke is defined by 
the equation 


24 ad h ad h COS iD 


and the requirement that /, shall always increase 
us qg; Varies By 12) and (10.2) 7, is then a singel 


valued function of /,, with derivative existing and 
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all / so that /, must be a s.v.d 
function ol By 1 nowevel q, is a s.v ad 
tion ol / so that 


Thus in a conditionally periody Staeckel svstem 


nonvanishing for 
Lure 
finally, qq, IS as. ad. function of 


anv coordinate Gg, IS a s.N d. function of he corre 


Spotl dl hy 


4. The »,’s as Functions of the Angle-Varia 
bles u 


evele of 
trip 
evele of a librational coordinate, we 


| J ina 


action and nngel variables 


If we now let nuh imcrease of T bye one 


rotatiol al coordinate and a sing round Irom 


ad to i bye Ole 
iv define thr 


by 
J pdg 14 


where 1 is now to bye considered i 
qs ict ott J’s rathet! than ol thre ( s and the as 
It is well know1 10, 11] that J, and w, are canoni 


cally conjugal so that 


- 


If Wie also define t] 


B » >. uf ‘ 
=F O-/; Oa “5 
where 
J. Oa 19 
creme / then le if to 
ip > / »() 
—s - 
But, by (7) ane iv 
] 
IB x O} , ~ O ia ms 9 
a= Oa 9 Il da i 
ol 
iB > oO qj | 2° 
r\ | t) na 4} Also. by 14 14 1) nie 4) 
pg q / 23 
lf we how introduce thre inatrix & thre niatrix 
q) P and tiv row” thatrices dv di ; : di 
and dw lu dw we find from (20) and 
> 
lrp=—duQ 24 


For a nonsingular Staeckel svstem m, > and & 
all eXist at every point ot (/, so that 
di Awd 25 
ol 
O Ou Qd eA 


Thus each exists and is single- 
valued evervwhere in @ 
by (14 Thus, by (7 
real, then WU 
nonreal 
olthe ws If all the w’s are real 
all thr ”s are 
tion Thus thre 
totality of all 
ps ure all real, also corresponds exactly to the sect 
lt therefore 
exists and 1s 

Thus each 


| Li uw 


derivative O Ou 
Now the J's ure all real, 
and (15 if the ps are all 


and the w’s are all real f some of the 


hen a Is nonre il ahd so are sole 
it then follows that 


should have 


p’s mre 
real, else we a contradic- 
corresponding to the 


which the 


domain ‘J, 
real values of the g’s for 
/ 
of all possible real values for all the w’s 
follows that derivative Ov,/Ou 
-Valued at an 
a s.v.d 


each 


single point il P-S pace 


nitst iy function ( In See 


tol } however, we showed that each g Is a s.v.d 


function of the corresponding Thus for a condi- 


tionally period nonsingular Staeckel svstem each 


isa s.v.d. funetion f,(w 


5. Periodic Properties of q,—/,(u , w,) 


By 7 ina Ld 


If now each coordinate q;, goes tt rough an integral 


of eveles, then by a familiar argument 


‘ a) ‘ 
Au y: " Op OJ lq S pag n 
— 7 i O./ — : i 
rs 
Thus if each dy. goes through exactly m, eveles, each 


Note however, 
alwavs 


Concel ned 


w, increases by the 
that simultaneous imecreases are not 
physically this 
only with the mathematical properties of the func- 


integer 
such 
possible section is thus 
tions 7,00 , W 

Sut we are really interested in the inverse problem 
eT om and 


Now the 


because of 


where each w has increased by an intes 


What has happened to the qs 
q's are uniquely determined by the w’s, 
the single-valued In the 
preceding paragraph where eacl 


we ask 


situation of the 
librational 
and each rotational 
angle variable 
Since the w’s determine the 
that whenever Aw 

each librational coordinate returns 


property 
coor- 
dinate returns to its initial value 
coordina ts gq, Increases by 2x” each 
u Increases by nN 
q's Unique lv, this has the result 


n j | n, 


to its initial value and each rotational coordinate 
gq, thereases by Dri 
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Thus, in the inverse problem, whenever we are | m’s. Note that the solutions for ? have no upper 


given Aw=—m, k=1, .. ., n, we find that each ro- | bound. 
tational coordinate g; must go through exactly m To apply this theorem, consider only those values 
cycles and that each librational coordinate g; must | of the time interval 7 such that »,7=/, where 7?’ 
go through some integral number of cycles, r,;, say. | is an integer that satisfies (36 In this time each 
By (28), however, we then find r;=m,. Thus when- | w, increases from its initial value w,(0) to a final _ 
ever the angle variables w, are all increased by | value given by = 
integer amounts Aw,=m,, each of the functions 
Ge=he(Wi, . . ., We) must go through exactly m, | w,( 7) =w,(0) + Pé,, 7 
cvcles. 
by (35 But by (36 v 
6. Mean Motions tne 
P: m Ne, n rf * 
If in a time interval 7 the number of complete | go that 
nnlis Neheied through by any coordinate gq is N,, w,(T)=w,(0)+m.+n 9 
the corresponding mean frequency n, is, by definition 
: As v T P takes on those large! and larger Integel 
Ne _lim N, T). (29) values corresponding to solutions of (36), each 
[> o approaches zero, by (38 Then, by (39), there 
exist infinitely many values of 7, with no uppe 
if the limit exists. We shall now prove that | bound, at which the orbit in w-space passes arbitrarily 
Ny=¥,=0a,/OJ,, k=1,..., n, for any conditionally | closely to points where Aw,=m,, k=1, . . ., mn, the 
periodic nonsingular Staeckel system. m,’s being solutions of (36 
To do so, note that if »,, . . ., », are all commen- If the initial g’s are given by 
surable, there exist a positive v) and positive in- 
tegers m,, . . ., m, such that qx(O) =f, [w, (0), .. ., wp (O)], c=, ,n), (40 
Vv; MpVo, (k Saas | 8 (30) then the values of the q's at anv of these times 7’ are 
given by 
where we may choose » to be the greatest common ( 
divisor of the »,’s. Then, from (16) and (30), during | q:(7 fw, (0 m ". , w,(0 m n I 
the actual motion, . 
pon] N 1] d 
w,—= w,(O) + myvol (k=1,. _n) (31 
As we let / P vy, ussume those larger and large! 
and in the time interval 7=1/%) we have values already referred to, the q’s then approach _ 
arbitrarily closely to the values W 
Aw,=m, (Pari, . . .7%). (32) " 
di | ‘4 fila, (0 my, , w,(0 mn" hy 
By section 5 each gq, goes through exactly m, eveles == ‘ 19 - 
in this time, so that in this case the motion is truly i" i 
periodic, with period 1/%. The mean frequency of | ,, oll ; ' tt] P 
qi is thus This conelusion fo owe Irom (3S) and the singk 7 
- m,/T'= mp a (32 valuedness and differentiability of the functions / a 
, : ; (‘omparison of (40) and (42) then shows that the ly 
If the frequencies », ..., ¥, are not all com- values qi(T) correspond to aw Mr, I l, ; 
mensurable. we mav let and are thus, by section 5, the values that would 
: reached after each qi had gone through exactly ” 
£,=p,/9, (ke Liewel (34 cycles. Now, by the definition (29), it follows that 


the mean frequency 


and then at least one of the é’s will be irrational. 


Then by (16) and (34), during the actual motion " lim (m, r), 3 
W, = Wx (O) + Exnit. (35) | if the limit exists. But m,/T=»,m,/P and, as 
T—o., lim(m,/P)=¢é,. bv (36). Thus, 
We now use a theorem of Dirichlet [12], which states 
that if the set of real numbers &, .. ., &, has at Np vee —=V_, 14 
least one irrational member, then the system of 
inequalities by (43) and (34). 
' | “Thus. for each coordinate qe of a gr gro 
&.—m,/P|\|<P a (k=1, ,n) (36) periodic nonsingular Staeckel system, the me: 


frequency n, is equal to the corresponding 7 
has an infinite number of integersolutions for Pandthe | mental frequency »,=0a,/0-J,. 
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7. Appendix 


Oo} 


dinates The 


corresponding Staeckel matrices and their inverses 


suite llite ol bits. 
oblate 


theories ol appropriate ( 


ure spherical Ol spheroidal 


ire, Wa sin & cos o, ¥ rsinégsind cos @ 


COs @O 


or.ita ( $ | 


easily written 


p 


hiost 


p 


Kor spherical coordinates @® bs 


down from the expressions for p nnd in the 
/ 


Kepler problem 1} with replacement Ol ab and a 
by 2a and 2a For oblate spheroidal coordinates 
® mav be found by cotnparing eqs 53) and (59.1) of 
} with eq (6) ot the present papet 

Note that ® or ® could fail to exist only when 
in 6=0 or when 7» | This could happen only 
when the satellite vroes over a pole and thus only in 


7 polar orbit Such a singularitv ina polar orbit 
however, is to be expected since nonsingularityv of 
n Stuer kel svsten leads to the q's being differentiable 
functions of the w’s and thus of the time In a 
polar orbit, on the other hand, the right ascension 
7 ® is a discontinuous tunction of time being 
constant except ut polat Crossings where if changes 


\ 
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Some Boundary Value Problems Involving Plasma 


Media 
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1. Introduction 

\ neutral mixture of positive and negative ions and electrons is often described as a 
plasma Kor example flames ruscous discharges strong shock Waves, and the ionosphere 


ire Various kinds of plasma. The behavior of electromagnetic waves within and in the vicinity 
of plasma is of great current interest 

Mainly because of the presence of the tree lectrons the plasma is a dielectri The col 
lisions of electrons with the tons cause dissipation of energy and thus the dielectric is lossy 
A radio wave propagating through the plasma is thus expected to be attenuated Furthermore 
In the presence of a steady magnet fie ld, the plasma becomes aMisOLropLle so that the dielectric 
eonstaunt is of tensor torm and thus 


rH reneral propagation is not reciprocal 
It is the purpose of the present paper to consider the interaction of electromagnetic waves 


wha plasma lor a special class ol two dimensional problems The reometry is chosen so the 
wave propagation is essentially transverse to the longitudinal axis. Such a restriction permits 
exact solutions of a number of relevant boundary value problems to be obtained Since radio 


waves, particularly microwaves, are being currently employed as a diagnostic tool in probing 
plasma, it uppears that results given below can find application 
An excellent introduction to the theory of propagation of electromagnet waves In plasma 


has been given recently by Whitmer [1] { more comprehensive treatment is found in a 
recent: monograph by Rateliffe [2] The dielectric behavior and the molecular properties ol 
it plasma have been dealt with by 2 number ol nuthors > 3] In parte ular spit hel 7|, has 


eiven a verv thorough discussion for fully tonized cases 


2. Dielectric Properties of a Plasma 


since the problen s to be discussed in the following sections deal with two dimensional 


reometry, it is desirable to take the axis of the Cartesian o1 the evlindrical coordinate system 
in the direction of the applied magnet field I In this case the dielectric displacement 
Dis related to the electric field E by the relation 
> > (1) 
DD € I 


whi re (e) is the tensor cielectri constant Kor an implied time frac tor exp wf), it has the form 
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The quantities ¢’ , e’’, and g are functions of the density of the electrons and the ions and the 
frequency of collisions between them. They also depend, of course, on the strength of the 
applied magnetic field //) (see appendix 
The case usually considered is where the electromagnetic forces only influence the ele 
trons. Furthermore, the motion of the ions is commonly neglected. For this situation, the 
properties of the plasma can be approximately described in a macroscopic sense in terms of the 
following quantities: 
Wy, the (electron plasma frequency 


wr, the (electron) gvro frequency 


v, the effective collision frequency (for electrons 


The elements of the dielectric tensor are then given explicitly by 


l } 
€ Ww I Ww 

q WrWy/@ 
t 

€ W7 I ‘Ww 

? 
€ Ww ™ 
l ) 
€ I Ww 


generalized to include the influence of heavy ions by simply adding 


a summation prefix to the ratios on the right-hand side of equations (3), (4), and (5 Then 
be employed This approach Is 


The preceding can be 


in each term, the appropriate value of w, v, and wr must 
valid for a weakly ionized medium and has been employed by Hines [8 
plasma consisting of a neutral mixture of electrons. one type of heavy Ons and t relatively 


In the ease of a 


large number of neutral molecules, the elements of the dielectric tensor have the form 


, . A 
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where 
P m mass of electron 
7m : 
m mass of ion 


and », and »,; are the collision frequencies of the electrons and the ions, respectively, with thy 
The effect of collisions between the electrons and the ions Is neglected 


uncharged molecules. 
the 


plasma and (electron) gyro frequencies. Since the charge of 


Wo and @r are the (electron 
7 
ure the 1On plasma 


ions Is equal and opposite to the electrons, it follows that Ze. and Uw 


and (ion) gyro frequencies, respectively 


It can be seen that if »,- w, the denominators in the second factors for e’ and q behave 


as 2w7—w*. Thus, at frequencies near the (ion) gyro frequency, the influence of the ions 
may be significant even though the mass ratio 4 is very small However, when collisions ar 


not negligible, the effect of the ions is usually quite small In particular, if: Ww and l 


this fact can be demonstrated by rewriting the elements of the dielectric tensor in the form 


( pw /w wrt) 
f] WW |W 4 > a 7 10 
€ W7 Vv; Tiger v 








Ao 
since the ratio W7 l “Wy 


: is of the order of unILY In most cases of practi al interest, 
The square bracket terms may be replaced by unity since rr is of the order of 10 


When the plasma becomes strongly ionized, the situation is more complicated. The 
formalism is given in the appendix where the relative number of ions and electrons is unrestricted 


However, in the case of a fully ionized plasma the elements of the dielectric tensor are given by 


9 10 
€ | m7 A mi Mb ] i 
1+ (12) 
€ A 
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| A (14) 
€ vlw—iv(l 7 
where 
\ nee p) Ww Ma aes | pL) 
As before, « is the ratio of the masses of the electron and the ton Also, wy) and wr are the 
electron plasma and (electron gyro frequencies, respec tively In this case, however, v is the 


collision frequency between the electrons and the tons which have equal and opposit« charge 


It e: 


n be immediately seen that if 2 is set equal to zero, these expressions reduce to the 


form of equations (3), (4), and (5), which are derived under somewhat different conditions It 
can also be seen that equations 12), (18), and (14) have the same form as (6), (7 and (8) if 

0 However, when collisions are non-negligible there is a fundamental difference 
between a we akly and a fully ionized plasma This Isa consequence of the coupling between 


the equations of motion for the electrons and the ions (see appendix 
Under the very reasonable approximation that 1+ can be replaced by 1, it readily follows 


n be written in the form 


that eqs (12), (15 and (14) ¢: 


€ L (WI au ~ 
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The quantity y could bye des ribed as an effec tive comple x collision frequency It should be 


noted that only the elements e’ and (j of the dielectric tensor involve - the remaining element 


«’ is not changed 





It is now evident that the motion of the ions, for a fully ionized vas, can be neglected only if 


a 


Ma Va 


Kon low lrequene ICs and Ol low collision frequencies this « ondition Thay well be violated. 
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3. Field Equations 


Maxwell's equations in a source free region with a (tensor) dielectric constant (€) and a 


(scalar) magnetic permeability uv are given by 


(ewok curl H (1S) 
and 
ball curl E { 19) 


It is desirable to write the first of these in the form 
iwh=(e~') curl H, 


where (e~') is the inverse of the dielectric tensor. It is not difficult to show that 


M ik O ) 
€(e °°) iK M 0 , (20) 
0) 0 € ad 
where 
M=—45 a md K=-Ge 


This formula is quite general. The only restriction is that the z-axis is to be taken in the 
direction of the impressed magnetic field. In the following it is assumed that the fields do 
not vary in the z-direction. In terms of Cartesian coordinates, Maxwell’s equations are 


then written 
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where N Mo/ M- 
It is a relatively simple matter to eliminate the transverse component of the fields from 


the preceding six equations. Thus 


eo, & , is 
latoetaan | die sale 
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Ou OY Ne a 

where I Eni us Jn wavelength These latter two equations ure valid only il the magnetic 
permeabilits and the elements of the diele« tric tensor are constant for at least a riven region. 

The fact that 77, and F individually satisfy a wave equation means that any solution to 
our problem can be regarded as the linear combination of two partial solutions. In the 
first of these, # 0 and in the second, // ) Thus without any subsequent loss of gx nerality, 
attention can be restricted to these cases It should be emphasized that this decomposition 
nto independent partial fields is valid only when the derivatives with respect: to ure zero 
As we shall see, the solution for 77,0 is relatively trivial since the constant magnetic field 77 
then has no influence (at least within the limits of mnacneto-lonic theory 


4. Reflection Coefficient for a Plane Boundary Between Free Space and Plasma 


With respect to the Cartesian coordinate Svstemn, a homogeneous plasma occupies the 
space y >0 and while y<_0 corresponds to free space, the constant and uniform magnetic field 
is parallel to the z-axis. A plane wave with harmonic time dependence (i.e., exp (iwt)) is 
incident from below as indicated in figure 1] 


The angle of incidence is @ (measured to the negative y-axis) and the wave is polarized 


such that Its thagnetic field has only a component ih the ~<lirection, denoted i Thus, 
Hl ho exp ikC'y) exp ikSx), 29 
where (¢ cos 6, S=sin @ and hp is a constant 


Since the refle ted field TT Isa solution of the free space Wave equation and is to have the 


sume dependence with z as the incident wave, it must be of the form 
HTs'=hy R exp Gk Cy) exp ik Sx). 30) 


where FP? is by definition the reflection coefficient 


With similar reasoning, the solution for the plasma (i.e., y>0O) must have the form 


HT,=f(y) exp (— uk Sa), 
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where f(y) is some function of y. Furthermore, since //, is to satisfy eq (27), it follows that f(y) 


must satisfy 


o* I 
, k2( — SS? f(y)=0 (32) 
dyt” VAIN ) bi 
Solutions are of the form 
exp [+ 2k [((.MN)~'—S?*)]'2y]. 3) 
Since the plasma extends to y= + ©, the solution corresponding to the negative sign in the 
exponent is pertinent. Therefore, the transmitted wave has the form 
H,=hy T exp |—tk [((MN)-!—S?]'2y] exp (—ik Sx), (34) 


where 7 is by definition a transmission coefficient. 
The boundary conditions are that the tangential components of the fields in the free 
space and in the plasma are to be continuous at y=0. Continuity of /7, leads to 


1+R 7. (30 


and continuity of £,, by virtue of eq (21), leads to 


C 1—R)=T [[((MN)-!—S?]'? M+ iKS]. (36 
Thus, 
R - . and 7 an , (37) 
C'+& C+ 
where 
A=M [(MN)"'—S?*]'?+iKS. 38 


For an electron plasma where the motions of the ions are neglected, it is possible to write 


where L (yr iw) w we, and Vv Wr W/ Wh We have also set N l (1.e., vi vi although a plasma 


A in the following form: 


may be slightly diamagnetic. 

The reflection coefficient, essentially in this form, was derived by Barber and Crombie [11] 
where the homogeneous electron plasma was to be an idealized representation for the ionosphere 
Because of the assumption of a purely transverse magnetic field 7/,, the horizontal direction 
of propagation is along the magnetic equator. For propagation from east-to-west, S is 
positive, while for propagation from west-to-east S is negative. y is then a positive real 
quantity. 

For applications at low and very low radio frequencies, : w so that to a good approxi- 
mation 

L~w/w where W Ww) 


Some numerical results based on eq (37) are available [12]. 


5. Reflection from a Stratified Plasma 


We shall now undertake to generalize the previous result to a plasma medium which is 
stratified in layers all parallel to the free-space interface at y=0. The situation is shown in 
figure 4 where Pr parallel lavers are indicated. The pth laver from the bottom is of thickness 
/, and its electrical properties are described by M,, N,, and K,. The index p ranges from 1 


to P. It should also be noted that /p= o. 
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Again taking the incident wave to be polarized with its magnetic vector parallel to the 


uXis, it is seen that the field for y< 0 has the form 


I1 hy [exp ke ('y I? exp ke C'y)) exp hk Sa +() 


The problem is to find an expression for 2 which involves the properties of the individual 
lavers It is possible ot course to formally extend thre results for the senu infinite case by 
writing the solution in each lave as a linear combination of the two elementary forms given 
by (33 The two unknown coefficients for each laver are then found from the two boundary 
conditions at each plane interface. The resulting 2.V linear equations may then be solved in 
a straightforward but a very tedious manner for any specified but finite value of The 


result init solution can be found ! i nore svstematic wav the analogy wit! ™ helk inoff’s 
13) theory of nonuniform trans lission lines is exploited We use this method her 
The wave impedances for the pth laver are defined by 


, I: ° I: tla 
A H and A I] 11h 


‘he superscript signifies that the fields varv with y according to the factor exp Boy 
where 
3 j VJ, N S 
wherens thre superscript signifies that thre fields Varv with y» according to the tactor exp 
8, Y n the present case, the superscript signifies a wave traveling in the positive » 
direction (14 iwav from the interface) and the signifies a wave traveling in the negative 


j direction 


Kron eqs tha ind tib Lt read lv follows that the wave inpedances mre 


A n, (VM, B k,S $2 


A n, (.V, B kK, S 13 


— 


Ln! € ~1?70 x The index p ranges fron 


Because ol the quantity kK 
it is seen that A and AK are not equal as thev would be u hn tsotropr mediun 


The reflection coeflicients at the intertace between the Q L)tt and thie Pt lavel ure 


now ck fine a by 
Hi] I: 


Ii i and / hI ’ $4 
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where the field components are evaluated within the (?—1)th layer. Thus, 


Kp_1—K} LK} L/W} 


a J. > i | ry , =, @ (45) 
Kp-itK>, ” l/K> +1/K; 


Rp 


The reflection coefficients at the interface between the (P 2)th and (P? L)th laver are then 


R K > Z; - 1/K; 1/Z, 1 
, > ; aii / ; { ») 
a K p-2+Zp . l/K> W/Z p 
where Zp_, is the impedance seen at the (?—1)th interface. From analogy with transmission 
line theory 
, ‘ 1+rp_; exp (—?72B8p_il = 
Le 1 A P s | . ‘ \ 4; } 
1+Rp_, exp (—i2Bp_1, 
where rp_, and Pp_, are given explicitly by eq (44 Now, in general, 
; -, |+r, exp | 28,1 
Z K, : ; (4S) 
; 1+ R, eXPp | 128, l 
so that the process may be continued until p=1, whence 
1—-Z 
R K 7 (49 
K,+2Z, 
Finally, for the bottom interface 
Kj —Z 
R ; (50) 
Ko +Z 
which may be rewritten 
C—& 
R 51 
Cra 
where A=Z,/n) since Kj =K No C. 
For the special case of a two-layered plasma (i.e., / ©), the explicit expression for A 


becomes 


> ( +28, 1,) 
A=—(M, 6.+iK, S) | 4 exp 2b 59 
1+ PR, exp (—? 2p; 4) 


with 
R (M,6,+7K,S)—(.M.8.+7K,S) a 
(Mb: —iK,S) + (M28,4 (KS ony 
and 
», — 1 /( MB +16 S) —1/(Mo8,+ 1K, S) - 
' 1/(\)8,—71K,S)+1/(M,8,+71K,S) 
The limiting case of a homogeneous plasma is recovered by letting /,-> ©, whence 


A=M, £,+: K;, S, 


which is identical to eq (38) after dropping the subscript 1 


6. Scattering from a Cylindrical Plasma Column 


In certain applications the plasma may be in the form of a cylindrical column. Examples 
are the ionizaticn associated with meteor trails in the upper atmosphere and the ionization 
associated with the shock wave emanating from an exploding wire. In this section expres- 
sions for the reflected or scattered fields are derived under the assumption that the ionized 
column is infinite in length First, the column is assumed to be homogeneous, but later the 
solution is generalized to allow for variation of the plasma properties in the radial direction 

Choosing a conventional cylindrical coordinate system (p, ¢, 2), a homogeneous plasma 


column occupies the space p<a. <A constant magnetic field, //,, exists through the plasma 
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column and it is directed along the z-axis. The medium surrounding the cylinder, p >a, ts 
taken to be free space The source of the field is taken to be a line source located at p py and 
d= and running parallel to the z-axis and the plasma column The line source radiates a 
evlindrical wave polarized so that its magnetic field has only a z-component 

The adoption of a line source rather than a plane wave source has the advantage of 
obtaining a more ceneral solution When p ©, the incident wave has a plane front but in 
most practical applications the incident wave front ts curved corresponding to a finite value 
ol p 


The primary or incident field is given by [14] 


ewl 


i] H'? (kp 55 


where / is the strength of the line source (actually, it is the magnetic current), and where 17 


is the Hankel function of order zero of the second kind, and 


f aD COs (@O @ 
En lovi Yr al idditional theorem for H kp . this can be written 
| 


I] 


for p< po, Where // and J... are Hankel and Bessel functions of order m. respective ly When 
p>p». kp. is to be interchanged with kp. The ¢-component of the primary magnetic field is 


t| hn viven by 


where the prime indicates a differentiation with respect to the argument kp 
Since the scattered or reflected field 77 outside the column is a solution of the wave equa- 


tion, it can be written in the form [14] 
/] 
where B,, is an undetermined coefficient The @ component of the scattered electric field is then 
hk >) BL? (kp) kl kp )e a0 

For the region p< a, the maenetic field /7, satisfies 


1 oO O 1 Oo j 
I dp” dp p? Op? M yl] 4 * 60) 


which is just the wave equation Therefore, we mav write 
J | 


enw / 


H > A I] kp J pie 61) 
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where 8=k(MN)-~"', and A,, is an undetermined coefficient. Since 


——s = oH, 
wh,—1k —M , (62) 
pen . poo Op 
it follows that 
E. H > AnH? (kos) | —™" J »(Bp)-4 MpI(30) | ™(o~<¢ (63) 
n=— p 


The boundary conditions at the surface of the plasma cylinder may be written 


H,—H'=+ H: , 
E kim KE he a. (64) 


Application of these leads readily to 


J, (ka) + B,, E75? (ka) 
J (Ba) 
and 
(*/) > Jn (Ba) mK J),(ka 
we N J,(Ba) ka J,,(ka) J (ka - 
. ( M ) 2 J), (Ba) m K HH?’ (ka I ka 
J ,,(Ba) ka HH)? (ka 


This is the exact solution of the problem. If the constant magnetic field //, is removed the 
results are identical to that of a dielectric cylinder in the presence of a line source [15]. In 


this limiting case, K=0, M—e/e and N= uo/u 


7. Scattering from a Cylindrically Stratified Plasma 


The preceding results are now generalized to a cylindrical column of ionization which 
consists of ? concentric layers. The situation is illustrated in figure 4. The incident evlin- 
drical wave again emanates from a line source at (po, d) and the field is to be observed at (p, ¢ 


The wave impedances may now be defined by 


; Ei. ; | 
K, —, and K : ; 67 
H H 
where +and—refer to the two independent wave solutions, proportional to J,,(8,) an 


HH (8,pe) in the pth laver. In view of the equation 


.-- OH 71 
/€ wk KK V/ : ’ ON 

poo Op 

it follows that 
K .Pmk, (M 2 J'.(B,p - 
me tp AN J m(Byp ne 
Pa 
FIGURE 4 1 line source in the pre 
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and 


: .PmkK, (M,\'? H?’'(8,p) - 
Kz,.=—i| F-(H) Ae gee? | (70) 


The index p ranges from 1 to P to signify the appropriate concentric region in the plasma. 
The reflection coefficients at the interface between the pth and the p— I th laver are then 


a. x ae 
 b..ara 


(71) 


R, 


and 

gE; Kt ...—iZ.. a 
Ejm I/K, |, aig 
where Z, »p IS the impedance at the interface between the pth and the p— I th layer Ln , 


mav be expressed in terms of the reflection coefficients 7?,,, and r, by again making use of 
. | } p 9: 


Schelkunoff’s nonuniform transmission line theory [13]. Thus, 


T.-K Fon, phen, AS p Cp +1)Xen,o\Ep+1y Gp) , (73) 


BE a,) as 


Xn e(G , a,) ; ’ (40) 


HH a 
H a 


(76) 


Hl a — 
x a ~a ° \é4) 


HH; (a 


The x's ure clearly transmission factors whit h describe the fractional change of wave as it 


il 
propagates from one interface to the other within a laver The numerator in eq (73) 1s thus a 
measure of the electric field at the pth interface taking into account the transmission through 
thre laver, reflection at the pP 1)th laver and transmission back again to the pth lavel The 
denominator of eq 73) corre sponds, ina similar manner, to the magnetic field 


The specific form of the transmission factors is 


\/,8.J,.(8,a mK ,.J,,(8,a a es 
X a ] a ’ is) 
\/ ee D,a mi a (pa a 
V78,H?'(B8,a mK,” (B,a,)/a ” 
X a a - ) 
Vo8.H Ba mk, H” (Ba a , 
x e >a 0 
I a ~ 
, J 3,0 
tie 
. /7. pa 
a a Si) 
tl Da 
Then starting with the reflection coefficients /?,,,p—, and r,,,p_,, given by 


as Ke as and } Ans Ji ~ S2) 


} oF | oe V/Ke , 1/K; , 
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we may obtain P,,,p_». and r,,,p_2 Which in turn enables us to obtain 72, and *,,.) The 
process is continued until the outermost interface is reached, whence 


ee. < to Ke .—WZ ” 
bm. 0 "K- =~ VA nie ; kK Y ’ » 
where 
: . Jd (kp 
k, of mf cp " 84) 
and 
. . Hl : kp) 
) ten ‘Ty (kp) n 85) 


If the scattered field 77% is written in the form 


Ti see > Be (kp vii (kp)e $s). ‘aie 

it follows that 
> > m A ka . 
B, K, He ka) i) 


It Thay be readily verified that 7,,.. for the special case of a homogeneous plasma column 
(i. e., Set d2=0) is identical to the square bracket term on the right-hand side of eq (66 

For purposes of computation it is convenient to locate the line source at a great distance 
from the cylindrical column in terms of wavelength. The Hankel functions of argument kp 


and kp, may then be replaced by the first term of their asymptotic expansions since kp and 
kp >>1. Thus, 


ine cowl 21 2 aa | 
a 4 ( rkp ; SS 
and 
, enol ya) 2 te ’ — | 
lit~ 4 —) ‘ p »> Boe ay] (kp é x9) 


This form of the solution would have been obtained directly if plane wave incidence was as- 
sumed at the outset. It should be noted that the factor preceding the summation in eq (89 
is just the value of //** evaluated at the center of the cylinder 


8. A Note on the Other Polarization 


The results derived in the preceding sections are valid when the magnetic vector of the 
incident wave is parallel to the cylinder. The derivation for the other polarization, namely, 
that when the electric vector is parallel to the evlinder, is relatively trivial since the constant 
magnetic field HZ) has no influence. The dielectric constant is now simply a scalar and is equal 
to e’’. The explicit results for E-parallel polarization may be obtained from those of //-parallel 
polarization by making the following transforms 


H.-F, 

E,——H 
iE, <naof 
E —+>—H, 
E—+—H, 
M+ N 
N-> M 


* 


Also, of course, K=0 and now M = (e/« 
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9. Appendix 
In this appendix a rather approximate theory of electric currents in a partially ionized 
ras Is YIVel The us 1S supposed to consist ol electrons ol Hiwass 7 and of ( hare ae and 
positive 1ons ol muss 7 and charge ( Since the vus rust be, tO a close approximation ele - 


trically neutral, the number of electrons and positive ions per unit volume is the same; each ts 


taken equal to .N There are also N, neutral atoms or molecules per unit volume with mass ./ 
indistinguishabl from the positive ons The electrons positive 1ons ana neutral atoms are 
regarded us three gases moving independently The interaction between the electron and the 
Oh Fuses Is supposed to by smoothed mto a continuous force In fact only those lorces are 


considered which result from macroscopic electromagnetic effects on charge and from friction- 
like action and reaction between the free charges ana the uncharged background Only mean 
velocities and forces are employed, and non-linearities are avoided by a perturbation treatment 

Furthe usSsuniptlions, adopted here, are that the velocity of the Tus us u whole is zero, 
and that the vradients of the electron and ton pressures are zero The removal of these latter 
restrictions would require that the problem be treated on the basis of magneto hvdrodvnamics 
‘ 16] Such uth approach has been v1Ven by spitzel 7| for il wholly ionized Cus and Cowling 
for a partially ionized gas 17] 


The mean velocities of the electron Cas, positive Ion Gus, and neutral Cus ure denoted by 
? Ui, and ‘ respectively sumce m m the electrons lose, on the average, a quantity of 


momentum equal to thei mean momentum m7 / at each collision with the charged ons If 


the mean time between SUCCESSIVE collisions Is 3 ‘ then the momentum lost by elec trons in 
collision with positive ions per unit volume is Nom,(¢,—¢,)%. Similarly, the momentum lost 


by electrons in collision with neutral atoms per unit volume is Nom,(v,—v?,)v, where vo! is the 
mean time between successive collisions of electrons with neutral atoms. Remembering that 
the momentum of the mass as a whole is zero, and since the mass of the electrons ts negligible, 


it follows that 
Noma,+N,m r, 0), (90) 
and thus ‘ at where a= No/N 
Now the electromagnetic forces acting on an electron are e(k- vi H where K is the 
ele tric field of the wave and H is the constant magnetic field superimposed on the svstem 


(It is assumed that //, is much greater than the Inagnetic field 77 associated with the wave 


The equation of motion for the electron gas may then be written 


Me d dt LL iwi ‘ kk Lal Il Hit itv - ar 9] } 


The equation of motion for the positive 1ons can be obtained in an equally simple fashion 


In this case, the momentum lost by the positive tons in collision with the electrons is 


\ Mi 


being equal and Opposite to the quantity uppearing above The positive ins, in collision with 


the neutral atoms of the same mass, lose half their momentum relative to the neutral vais 


Specifically, the momentum lost ts Vom av aiv, where v is the mean time between suc- 
cessive collisions of positive tons with neutral atoms Therefore, the desired equation Is 
madv,/dt Mi, tae My kk Mat H] Wi Voll my vy l a)/2 (92 


Equations (91) and (92) may be solved for 7, and v; in terms of # and other known quan- 
tities. The current density resulting from motion of the electrons and the ions is then equal 
to Noel . Noting that the displacement current is je@/? and the total current is ((e)ol, 


it is seen that the dielectric tensor (e) can be obtained from the relation 
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i(e)wHh =e, wE Nee(e, r,). (93) 


> 
Choosing H, to be directed along the z-axis, the dielectric tensor is found to have the form 


e’ iq 0 
(€) iq e 6+ (94) 
0 i ss 


The quantities ¢’, «’’, and g are functions of m,, m,, »., »;, v, No, Ng, and Hy. Explicit results 
for certain special cases are given in the body of the paper. To simplify the notation there, 
the results are expressed in terms of the positive real quantities w) and wy which are defined by 


we = N,e*/em, (95 
and 
Wr pol Te m, (96) 


Also, in discussing certain special cases in the text, the subscript on » is dropped. 

1 would like to thank Mrs. Eileen Brackett for preparation of the manuscript and Mr. 
Kenneth P. Spies for his careful reading of the paper. The present work was performed under 
contract CSO and A58—40 with the Electronics Research Directorate, U.S. Air Force Cambridge 
Research Laboratories. 
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A New Decomposition Formula in the Theory 
of Elasticity ' 


J. H. Bramble’ and L. E. Payne 


February 17, 1961 


In this paper a new representation formula for the spherical components of displacement 
in the theory of elasticity is obtained. This formula may be used to reduce mixed boundary 


value problems for the elastic sphere to standard problems in potential theory 


1. Introduction 


In treating boundary value problems in the theory of elasticity it is frequently convenient 
to represent the components of the displacement vector in terms of harmonic functions, and 
thus to reduce the elasticity problem to one in potential theory. Several such decompositions 
in terms of harmonic functions are known, [7, p. 342-343],° [1, p. 231], 
perhaps, being the Papkovitch-Neuber decomposition see for instance [7, p 331] 

Half-space problems have been quite thoroughly treated from this point of view, as in 


fact have the first and second boundary value prol lems for the sphere However, In treating 


[4], the most widely used, 


mixed boundary value problems for the sphere, the standard decompositions are not suitable. 
Such mixed problems could be simplified considerably if it were possible to obtain convenient 
representations for the spherical components ol displacement in terms of harmonic functions 
We obtain such a representation in this paper and demonstrate how it yields in a simple man- 
ner solutions to standard mixed boundary value problems for the sphere 

In this paper whenever an index is repeated in a single term, summation from 1 to 3 will 
be inplied Use of the svmbol f,, to denote Of/Or;, where sz; 1s one of the rectangular coordinates 
J Be. Ba. will be made throughout 

Let D be i three dimensional region which contams an open subset ol a spherical surface 
G. ©) The region D Is 


of radius a with center at the origin ol a spherical coordinate system (7", ¢ 


assumed to be TT portion ot an isotropic homogeneous elastic medium 
We suppose that )) does not contain the origin. Furthermore, we assume that every ray 
from the origin which intersects J), intersects it in a single line segment 


In J) the component of the displacement vector satisfies the equilibrium equations 
u au 0) 1.1) 


1—2c)~':where o is Poisson’s ratio). As far as our 


where a involves the elastic constants (a 
representation formulas are concerned we need put no restrictions on the constant @ other than 


that it be different from —1. However, all physically interesting values of a lie in the interval 
| 
> a x 
> 
I ‘ pt lin part by t U.S. Air I 1 Air | ort Ss itific Researcl Air Res Ly 
I nt ¢ ul | ( tract, Ne A F 49(638)-228 
University of Maryland, College Park, M 
Figur t ke tl teratu t! pay 








In this paper we obtain the following representation for the spherical components of dis 


placement in D: 


, G\) 
ryu a? 2G + (2—a)rdG/Ov a ° ( ° ) t/ ov 1.2) 
O/ O/ Oo 


P . 7) H] 
€ fp (Uy ed. (2+ 8a)e,,2,G a €, 22,4 —( H+, : ) 
O/ O07 


\ 
oH. d/ oll all | 
+ {: lr ‘= +27 >) (, “ | -} [ > | } Hé, rv, 3 


where ¥, G, and H/ are harmonic functions in ). Consequently all terms in braces in (1.2) and 
(1.3) are also harmonic. In (i.3) €; is the permutation svmbol 
| 0 if any two indices are equal 
€; l if is an even permutation oi233 (1.4) 
1 af ¢, is an odd permutation of 1 23 
2. Representation Formulas 
We first rewrite (1.1) as 
: I 0) BI 
Ea Carn Ur mT iC iu (Z.1) 
ou 
For a # —1 we see from (1.1) that w,,;is harmonic. It has been shown that the representation 
0 : OG 
Myy=(1+r 5, )(6G+4r 3° ), (2.2 
Or Or 
where G@ is harmonic, is valid in D) (see for instance Bramble [2]). Equation (2.1) then may be 
integrated to rive 
' 0G 
€: in Uy 3 (a+1)e« x,( 6G+-4r ‘) +X, (2:3) 
; or /,) 
where x is a harmonic function. We express x in J) as 
oO 
d+r—. (2.4) 
O/ 


Thus u; is the component of a vector whose divergence is represented by (2.2) and whose curl 
is expressible by (2.3). The component wu; may clearly be expressed as the sum of a particular 


solution of (2.2) and (2.3) and the gradient of a harmonic function, i.e 


’ OG a. , ) i 
Uj 22] G+(a+2)r |_rl 21 t+-Da) G ar 9G. Hear, ty (2.5) 


where ¥ is harmonic in DJ, 
Using (2.5) we now form the expressions for zu, and €;,2;u, which are respectively the 


scalar product and the vector product of the radius vector and the displacement vector. 


. ' 0G OG 
I Uy; / 26 +-(2 a)T —ar?() ) 1, (2.6) 
Or O/ O/ Or 


2 9.1.2 ’ O ' OP — 
€(j.0gh =r" —(2-+-32) €:).35G,, ar> €15 0%, ?P rir > mL ILA (2.4) 
/ / 


and 


<< 


Equation (2.6) is precisely the desired expression (1.2 In order to obtain (1.3) we substitute 
| | | 


lol ? in dD see Bramble 2] 
e=H+2r! 2.8) 
O 
where H7 is harmonic We then write 


On o// o/ oly) ofl ( oH am 
On [4 Fp ar 2 (OM) S (3) ees) - 


As noted previousl\ the term in square brackets tis harmon Equations ey 2.8) and (2.9 


vield immediately the desired representation Lo 


In obtaining oul decom oOsitions Wwe have free uently made use of the fact that in the 
| | 


: O O 
region inh question any harmonic tunetion could be re presented aus ( l+ay > )( |l+a) > ) g 
/ 


where is harmon We have used only representations in which a and a. were non negative 
sSramble [2] has shown that such a representation is valid in D regardless of the sign of a, and 
a However, for a, and a. non-negative it Is easily verified that such a representation 1s also 
valid mb anv region whose bounderv ts star shaped with respect to the origin ( onsequently 


1.2) and (1.3) are valid in any such star shaped region In the next section we use the repre- 
sentations (1.2) and (1.3) to solve certain boundary value problems for the elastic sphere 
If the elastic region is exterior to a surface C which is stat shaped with respect to the 
origin and if the components ol the displacement vector are OU us / it can again be 
shown that the decompositions 1.2) and (1.3) are valid therein 


3. Mixed Boundary Value Problems for the Sphere 


Let vw, Sutisiy 1.1) on the interior ot an elastic sphere On the surface a let 
a rid Gotu hG.¢ 
5.1 
L h € J j h é j / qd ed rea 
where 7 is the stress tensor, u is theshear modulus and where f and q are prescribed boundary 
data Equations (3.1) may be written equivalently as 
O oO 
( Y ] rou Aa\\a Gdo\£,u a Y € J € ry u AU 
O Ou 
@) 
h é rou h oh é u h é j r,u GAGwo 4 2 
0 , 
on? (i Using 1.2) and 1.3) we can actually obtain the following representation [ol the 


left hand side of (3.1) in terms of harmonic functions 


; + >) + } 7 
hn Ose I olf | o 2 (, oll ] 2 (, OH . 
(bur 2 +b. 2b )q aif ME par 2 (eM) Le (VE) VG 








From (3.3) we obtain expressions of the form 


A (> < ) @+4,(r 2) @+Ay C+ AG+Ay S (7 ©) 4 Ar =P), 


or 
een? 2{s (r 2 ) G+ By ~ 1 BG+ By oY \ BS (3.4) 
+{ «(2 S41) (7 2) (Br S24 Bat) —2 | r 2 (Br 4 4B.) | } 


[lr d) (mes ma)) 1 


valid throughout )). The constants A, and B, are determined immediately from (3.3) an 
P(f) and P(g,) denote Poisson integrals. 
Forming z,/’(g,) we obtain the expression 


" 0 \” e* . o\? oll z,P(¢ _ 
Be {(r2) H+(r 5) u} + B; {( =) H+r— }- ae (3.5) 


Under suitable conditions on B, and B, this expression Day be integrated to give an expression 
for 7. With known throughout D, equations (3.4) give two identities in the two unknowns 
Gandy. The function ¥ can be eliminated by taking the proper combination of r derivatives 
of the first of equations (3.4) and r derivatives of the expression formed by operating on the 
second of equations (3.4) with « miki > - This results in an equation involving a_ linear 
Lea 

combination of r derivatives of G. If this equation can be solved for G then the insertion of 
G into the first of equations (3.4) vields the expression for ¥ which completely determines the 
displacement vector in D). It is a simple matter to obtain a formal solution for G, for we have 
an equation of the form 

: O 

: ( )@ F. 3.6 


r 
i=0 or 


We expand the known harmonic function F in terms of spherical harmonics, i.e 


N ra) i ’ 

=I . » Dm = 
Sa (r ) G . S* r"P™ (cos 6)[A,,, cos mo+ B,,, sin mod 3.4 
i=(0 or n=0 m=0 


The solution to this equation is 


x n N l 
G > > 2a a r"P™ (cos 6) A,, cos Mmo~ B sin mo 3.) 
n=) m=0 i=0 
provided 
N 
S*an' +0 (3.9) 
i=0 
for n=0, 1, 2, . . ., or any positive integer. For a given a, (3.9) imposes certain restrictions 


on the constants a,, do, b;, and by. 
As a particular example consider the case in which a, and ), in (3.1) are zero. The constants 
a, and 6, can without loss be taken as unity. Thus equation (3.4) becomes 


(r2-1)y P(f) (3.10) 


, G 
a'| - 2a(r°-) G+2(4a-+ yrs + 2(3a 1G }+2 ° or 


, i 
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mati 
) 0G : oH 
O° €; pI —al (2+-3a)G €, xt a*( H ; ) 
Ou 07 Oo 
oll o\? oll 
+ +2 H Pi 3.11) 
{ «|: O7 (r5") ] A(r>) d 
where P q and P /) denote Poisson’s integrals If we take the scalal product of the radius 
vector and (3.11) we obtain the differential equation for I] 
2) oH r,P(q;) 
} i] / - : et 3.12) 
ey) Or _ r*—a‘ 
which yields the result 
p Pi 
i gs ett da 5.13 
. —) a 
It is possible to »dd an arbitrary constant to the mght hand side of (3.13) but since only deriva- 
tives of //7 occur in (1.3) we may take the constant equal to zero. We now operate on (3.11) 
O 
with the operatol é j to obtain the expression 
Ou 
O . 0G @) 
r— < rl ~¥—(2+3a)a°%G—a’ar — €ssa2 P(q ; 14) 
O/ O7 Ou : 
which vields for ¥ 
a : Oot op , 
¥ a E sa)G r/ | / | p / pje r ou P(g \dp- » ) 
Again we have set the arbitrary additive constant equal to zero in (3.15) since only derivatives 
of w occur in (1.2) and (1.3 Substitution of (3.15) into (3.10) yields the result 
O re) : ' 
(, >) ( )() >) ( )G ta-(a I 516 
where 
( ) a ‘ 22a 23a") ’ } a l 
3.17 
( I—a 17 —22a—23a*] tia l 
and 
oO Se oO 
R P / 2(7 | )| ‘| € J Pla do | s.1S 
O07 . or : 
For a >1/3 (which corresponds to the range of ¢ which is of physical interest, 1.e., l< o< 1/2) 
c, and ¢ have real part less than unity. Thus if we operate on 3.16 by rd/Or we obtain 
) o od , OR 
(, —¢ )(y ri )y ; ta*(a+-1) r : 3.19) 
O07 O7 Ov or 
an expression which can easily be integrated for G, giving 
(0.0.0) ? *% *? OR 
G 4a*(a+1) { : T p ‘| p ‘| p, ' : dpodpidp . 3.20) 
ee Ji , Jo ——- OP» 





If the indicated integration is performed we obtain for ¢,#c, the expression 


R(r,6, : ’ e } 
G=|4a*(a-4 1)] . = ) 4 (6, C>) ; p of dp- ; | p of dp > (3.2 
C)Co ~ c Op ¢ Op J 


If e;=c, (which means that a@ satisfies the equation 17—22a—28e=0) (3.20) vields upon 


’ 9.3 is oR j OR 
G=|4a?c?(a+1)]7! qr pp“ log (r/p) dp—) p dp+R ?P- 3.22) 
; e Op Op 


The function y is now determined by insertion of (3.21) or (3.22) in (3.15), and the problem is 


integration 


solved. 
4. Concluding Remarks 


In a recent paper [3], the authors indicated how to obtain an explicit continuation across 
a spherical surface ( on which f and g,; of (3.1) were zero. The decomposition (1.2), (1.3) could 
also be used to obtain this continuation formula. It is also possible to obtain in a straight- 
forward way from (1.2) and (1.3) sphere theorems with various combinations of surface trac- 
tions and displacements vanishing on the surface of the sphere 

If in (1.2) and (1.3) @ is allowed to tend to infinity while at the same time aG@ tends to a 
harmonic function F’, then these decompositions may be interpreted as representation formulas 
for the spherical components of the velocity vector in the slow motion of viscous fluids See 


for instance [6, p. 522].) 
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Pointwise Bounds in the Cauchy Problem of 
Elastic Plates 
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et a elastic plate ccup i reg D wit ! lary ¢ () ij t > (” the 

fle« lope moment mad ir are meas ed l terms i thes data cr le 

s rt maxim flect i moment, | Vise | is for the va the 

lect and its cle " ‘ D ! tained Dhese ‘ | ire such t t the 

lata Sa i the load yl ! j air ft t pia the plate te it el t lec 

; j ; at « ! ; D wiy ' } 1 bta ed ; 4 

! ra j 
1. Introduction 

The following plate problen is considered Let an elasth plat occupying a region DD 
a thre ry plane he sub on ted ton load PILY normal to the plane ol the plat The deflec tion 


Sutishes the equation 


O O'v ?) 
Ar =( 2 )=p/d Ll 
du Ory OV 
where J) denotes the plate rigidity (see [3, p. 488 We assume that on a portion = of the 


boundary ol DD we are abl to measure the de fle tion, slope the sum of the two prim iple no 


ments, and the normal component of shear with known error: o1 equivalently we assume 


/ a I | d y 
. . 19) 
A / , A G,)(A G,)d ay | 
. . r 
where | I q ii al G ure thi measured values ol A and Ar, respectively In La 
the comma denotes differentiation and a repeated subscript i implies summation over 7 (7 
| 2 The A svinbol is used throughout Lo denote the Laplace operator 
In (1.2) then the a, are assumed to be small known constants We wish to obtain upper 
and lower bounds for the solution ¢ at any point in 2? 


It is well known that as it stands the Cau hy problem for the biharmont equation 1s 
not well posed see Hadamard [1 The solution does not depend In a continuous Way on 
the boundary data We leave aside the question of existence of solution and merely remark 
that in most physical situations in which one would be interested in considering such a prob- 
lem the question does not arise Kor instance, we might wish to treat an ordinary plate 
boundary value problem, but find that part of the boundary is inaccessible for measurement 
ol boundary data whereas on another portion of tI 
ov, Av, and O (A 
On oO? 


If inh addition to the Cauchy data we have uniform bounds for both and A throughout 


e boundarv it Is possible to measure ?P. 


D) then the solution to the plate problem will be seen to be stable This is in fact an im- 


mediate consequence of the results of Laurentiev [2 
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In order to simplify the procedure somewhat we decompose the solution v into 


v=0,+0’, (1.3) 


where v, is any particular solution of (1.1) defined in D and v’ satisfies 
A?v’ =0 (1.4) 


in D. We approximate v’ by a biharmonic function ¢ and set w=v’—g. Thus w satisfies 
in D 


A*w=—0. (1.5) 
Also bounds for the Cauchy expressions 


MN ude, Ne w,w ,de 


. J = 


€ (Aw)*do, € Aw ,Aw ,doe 


are known. In fact 


m1 [(v’ r0,—] )—(oT+?, f)\*do < 2a; 2| oT, I] \*do (1.7) 


with similar inequalities for 7, «, and e;. It is clear that if a; is small and @ can be chosen to 
approximate f-v, in mean square on Y then », will be small. A similar statement holds for mp, 
€,, and ¢. We shall show in this paper that it is possible to obtain an a priori bound for 


w of the following type at any point in P?: 
w( P)\?< Ay Mt~**{ ay + aon, }®! + ApM}~ 2! bye, + doe. }*, (1.8) 
where A,, A», a), az, 6), b:, are determined constants, 6, and 6, are constants satisfying 
0< 6.<6,< 1, 1.9) 


and VM, and M, are the uniform bounds for w, and Awin J). It is clear that knowledge of uniform 
bounds for v’ and Av’ implies knowledge of uniform bounds for w and Aw since w-r’ is a known 
function. 

With an inequality of type (1.8) it is clear from (1.7) and the corresponding inequalities 
for m2, €;, and ¢, that if the a; are sufficiently small and the Rayleigh-Ritz technique is used to 
make the integral expression on the extreme right small then we obtain close upper and lower 
bounds for v’(P?). This paper will be concerned then with the establishing of (1.8). 


2. A priori Bounds 


We consider a two-dimensional region D) bounded by a closed surface C. On a portion 
> of C, Cauchy data are prescribed. We denote the remainder of the surface (C-2) by ¥. 
(As will become apparent later the portion > need not a tually be prescribed.) We now choose 
the origin of a polar coordinate system at a point ?y) which has the property that there exists a 
circle K(r) of radius r and center at > such that the intersection of C and K(r) contains no 
point of =. (If two or more disjoint regions are formed by the intersection of K(r) and D, it is 
sufficient that the boundary of any one of these regions contains no points of 2.) We denote 
by D, the intersection of K(r) and D (the portion whose boundary contains no points of & if 
two or more disjoint regions are formed) and require that /’y lie outside of D,. We denote by 
R the radius of the largest circle satisfying the above conditions (for a fixed 5) and by ro the 
radius of the smallest such circle. As indicated in the previous paper [5] we locate /) in such a 
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way that A(7») is tangent to )) at some point of =. Finally we denote by w(r) the portion of 
the surface of A(r) which lies interior to ¢ 


As in [5] we introduce for any function uw the notation 


> J 
€) “u-dod 


f(y +Qy 


where do denotes an element of surface area on 2, and let 


Asoo [J[(32) +35) ate 22) 


2.1) 


vs 
In [5] we proved that if w satisfies 
Au=0 u<M in D, 2.3) 
then the functional 
rur, u loz | | p A (u)dp-+ K€ mee +4 (log 7 . (2.4) 
is a convex functional of log r in ro< r<— R, where 
; 3 > 1 De QF 
( -Al y2)7 I K = 9 h oe log (Ry) ) 27+ fi j (2.0) 


It was emphasized in [5] that the constants in [2.5] are usually not the optimal constants, which 
In anv spec i fic problem depend on the geometry of = and 1) as well as the location of P? 


An evaluation of the integral in (2.4 leads immediately to the inequality 


P | 
|p 'A(uldo—, wds | <r "etry log (R/ro)e (2.6) 


- 
a 


ead to the inequality 


} u*ds ry € R?, é 2BM » (qe, + Ke)? 2.4 
The constant B Is easily comp ited and in this case 
6=log (R/r)/log (R/) 2.8) 
Try the plate problem under consideration (1.5), (1.6) we now set 
Aw=u 2.9 
Then (2.7 is valid with wa replaced by Aw since we shall need it later we obtain now by 
integration the following inequalits 
ey , | l ; = 
Aw dardy Ts } ; r ars re) hi , é 
D é ~~ 
/ A; € ’ A»€ 
} V/ Ki€ Ko€ I ( VE ) Ji Rii 
pr? : 
2B 2 10 


Note that for 








which must hold if the estimates are to be useful, we obtain the cruder but simpler inequality 


- ; Fil 
(Aw)*dsdy+; (7?— r5)~ roe Ts (r?—72) R8ro en << BPM K€ Ko€ 2.12 
?- 

We wish to find a bound for u*ds. But once such a bound is known then pointwise 


bounds for w at points in D, follow at once; for we need merely introduce the biharmonic 
Green’s function r(P,Q) which vanishes with its Laplacian on the circle of radius r and centet 


at >, and has a fundamental singularity at a point ? in J). We then make use of the identity 
: a or ) = or oO 
w(P)=— w — (AT )ds— Aw <— ds— w— (Ar )—~ 4r+4w —- Aw)T |de, 
ot ots on J on J: on On Oo On % 
912 
where Y’ denotes the portion of © which lies interior to A(y (n application of Schwarz’s 


inequality then vields the desired result. 


In order to obtain the bound for u*ds we introduce the following functional 
U'(rw) lo { p | As(w) +. wawedarly | dp 1] n> lens ae | Nat ion} OF oe ) 14) 
eo Ta “Dp 
where 
1 unde, "2 ww, de. 
™ ; 2.15) 
"3 (Aw)*dady, ns=2 BR, M?P°-® (nye + K€ 
De 
for 
r<or< hick; B=log (R/R) log (R/r 2.16) 
and 
p>|, qg>p—l. 2.17) 


We demonstrate that for p and q satisfying (2.17) and the constants /;, /», /s, 44, and ©; suitably 
chosen, U’ is a convex function of r 
To this end we introduce the notation 


U(r, w)=log E(r, w)+- Cyr-*4 2.18) 


(where £ is the term in braces in (2.14)), and form 


ger a2 b{ (q+1—p)r**! *| A, w+ | | wud ey | f- 724 ( | li grad w * 
D 
+wAw|ds )} —p2atl—p [ 4. w) + wawedely | +29 2.19) 
”. 


As in [5] we obtain in a straightforward manner the inequality 


ee) NoNe, 2.20) 


a 
Min, + Mone < E—= rr 


°° 
w'ds— 5 (p—1) p w*dxdy—l3n l 
A 
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where 


) 
ii j 9 > , i 9 p . | 
2 >] } 


The constants / ind / are to be so chosen that 7 nnd w ire positive. 


$y Green’s theorem 


side | ‘ ae . —_ 
A.(w)- | | vAwdardy ” da DY de » 29) 


Thus by apply g Schwarz’s inequality in (2.22) and using the fact that the arithmetic mean 


jo 7 
‘ 


Ov l ” [ b <i Y 


Is not less than the reometlric mnean we obtain from 2.19) using (2.21 and (2.9 


ree p—I 
4 2 qr { - / ds + u = | uw lady lon lang Win Mon } 


where a, 6), b:, 63, and 6, are arbitrary positive constants. From the Rellich identity [6, 7 


it follows that 


‘ Ow Ow *, ow Or Ow Ow 7 or 
' {( O )-{ o ) ba . 3 Awdadyt2 "Oa Oa O do ae grad w de 2.24) 
~ } j x 0 Oo Ni : Oli 


Thus, for any 6,0, 
Ro) ee 
If ov Or is to be non-negative, then the coefficient of the term A,(w , u’-ds must be 
non-negative After insertion of (2.25) in (2.23) we find that this eneflicions is positive if 
q—1 l 5 >0 2.26) 








The coefficients of all other terms may be made positive by choosing () sufficiently large 
provided 
p>}. 2.24 


But (2.26) and (2.27) are just the imposed conditions (2.17). Hence for such a choice of the 
constants, U is a convex function of r~*. 
We now let 
6*=[(Ri/yo) ~*—1)/[ (i/o) ~*— 1). (2.28) 
It follows then that 


wp 20, a “88 _— =e ry - 
ec" E(rw)<[e" E(R,w)]'-* [eo E(rw)|*. (2.29) 


Using (2.20) we obtain the inequality 


1 . = (p—1) ; 
“LA "? w'ds+ PY wir ‘er! dxdy + My My + Mono + ls ngtl yn, 
<O,Mi-* {myn t msg.) +C2ME-© {eye t+ e2}*", (2.30) 
where the constants (, and (, are determined as indicated in [5]. We have used (2.10) or 
(2.12) for ns and the notation 
6’ =66*. 2 33) 
By use of Schwarz’s inequality and inequalities (2.7), (2.9), and (2.30), equation (2.13) then 


yields the desired pointwise bounds for w(/) at any point P in D) 

In order to obtain the desired bounds in ))— D, we first form from (2.13) w,,, Aw and Aw, 
For P a finite distance from the boundary of D, the differentiation may be carried out under 
the integral sign on the right and all of the resulting integrals converge. Thus using (2.7 
(2.9), (2.13), (2.29) we obtain the pointwise bounds in DP), not only for w(), but also for the 


quantities |grad w|, Aw and grad Aw (in fact estimates for any derivative of w If we now 
integrate these inequalities over a portion of a spherical surface w(r;), where ry<r and each 


point on w(7,) is a finite distance from 2, then we may regard w(r;) as a new Cauchy surface 
and proceed to define a new region J),,. continuing as indicated in [5]. We are thus able to 
obtain pointwise bounds for w and its derivatives at any point on the interior of D) 


3. Analogous Viscous Flow Problems 


Instead of interpreting our biharmonic Cauchy problem as a plate problem we could also 
have considered v in (1) to be the stream function in the slow flow of an incompressible viscous 
fluid [4, p. 541] (here p=0). However, in this problem the actual value of v is of no physical 
significance. The velocity components are expressed in terms of derivatives of v, and the 
physical quantities which we measure on = are the velocity components, the vorticity and the 
pressure gradient. In terms of v the vorticity is the Laplacian of ¢ and the pressure gradient 
involves the derivatives of Av. 

Hence we are interested in estimating only the derivatives of ¢ and not ¢ itself. At the 


same time we would like for these bounds to involve only m, €, and &, and not 7 In order to 


accomplish this we consider instead of (2.13) the identity 


w(P)—w(P) w w(P)| (AT )ds— ao = ds 
J wir on J a on 
i 5, 0 ow or oO 
si fw w(P) (4r)— AP+ Aw (Aw) 0 $do 5 1) 
J on on On on 


where P is the point on Y at which A(7») is tangent to ¢ 
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For points away from the boundary of D), we differentiate (3.1) carrying out the differenti- 
An application of Schwarz’s inequality \ ields the same types of terms 


ation of T on the right 


as those encountered previously except for the terms | w uw Ids and Ww w(P) Pde 


Using the arguments emploved in deriving (2.29) with w replaced by w—w(/’), we may bound 


the integral over w(r) in terms of no, &, « and Ww w P "doa Let us now define A by 
"lO 
(2%) do 
P Oa 
A=min - 5.2) 
among continuously differentiable functions y) which vanish at the origin The number A ts 


the lowest eigenvalue for a vibrating string fixed at the origin and free at l. Thus 


\ WT 1/ 3.0 


We now let l be the shortest distance on > trom PP? to either of the points ol ntersection of 


and S. and denote by ZL the length of Thus 


l< L/2 3.4 


and since wu w(P?) is clearly an admissible function WY im (2.32 the origin being taken at 
a Wi obtat using > and ,.4 
_ Ou - 


It is thus possible to obtain in J), the following bounds for the velocitv components w using 
| 


27 2.9 2.50 Rs and (2.1 

pu (CM ) (MM K€ Ko€ 6 
where ¢ and ¢ are constants and all other terms are the same as those occurring on the nght 
In (2.50 Note that as in [5] the constant V/ which is a bound for w w( P mav be replaced 


hy 1iconstant L Wile hi involves : bound for thre Kinetic em rey inh DD 
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